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1  Executive  Summary 

The  principal  goal  of  this  three- year  research  effort  was  to  enhance  the  research  base  which 
would  support  efforts  to  systematically  control,  or  take  advantage  of,  dominant  nonlinear 
or  distributed  parameter  effects  in  the  evolution  of  complex  dynamical  systems.  Such  an 
enhancement  is  intended  to  support  the  development  of  flight  controllers  for  increasing  the 
high  angle  of  attack  or  high  agility  capabilities  of  existing  and  future  generations  of  aircraft 
and  missiles. 

The  principal  investigating  team  has  succeeded  in  the  development  of  a  systematic 
methodology  for  designing  feedback  control  laws  solving  the  problems  of  asymptotic  tracking 
and  disturbance  rejection  for  nonlinear  systems  with  unknown,  or  uncertain,  real  parameters. 
The  main  technical  issue  which  was  overcome  was  the  lack  of  a  traditional  observer  -  or  state 
estimator  -  design  for  uncertain  systems,  caused  by  the  fact  that  traditional  observations 
cannot  detect  the  values  of  the  unknown  real  parametric  uncertainties.  This  limitation  was 
overcome  by  the  development  of  a  new  design  methodology  which  is  based  on  the  concept  of 
immersing  one  system  into  another,  a  topic  which  is  discussed  more  fully  in  this  report  and 
in  a  recent  book  by  the  principal  investigating  team.  The  idea  developed  was,  intuitively, 
based  on  an  attempt  to  propagate  just  the  information  needed  in  a  feedback  law  rather  than 
the  entire  state  and  parameter  values.  When  applied  to  the  Kalman  filter,  this  observation 
also  led  to  the  solution  of  a  fundamental  problem  in  stochastic  systems  theory,  with  applica¬ 
tions  to  speech  synthesis  and  voice  recognition.  A  patent  application  for  a  related  invention 
is  pending. 

Another  successful  research  project  was  the  development  of  a  systematic  feedback  design 
theory  for  solving  the  problems  of  asymptotic  tracking  and  disturbance  rejection  for  linear 
distributed  parameter  systems.  The  technical  details  which  needed  to  be  overcome  are 
discussed  more  fully  in  this  final  report. 
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1.1  Robust  Nonlinear  Control 

As  a  starting  point,  we  recall  that  one  of  the  classical  problems  which  involves  shaping  the 
response  of  a  system  is  output  regulation,  in  which  the  objective  is  that  of  controlling  a 
plant  in  order  to  have  its  output  track  (or  reject)  exogenous  commands  or  disturbances. 
The  nonlinear  regulator  theory  developed  in  [12]  is  valid  in  case  a  model  for  an  “exosystem” 
which  generates  the  exogenous  signal  to  be  tracked  and  the  exogenous  disturbance  to  be 
rejected  is  known  and  satisfies  certain  conditions.  These  conditions  currently  prohibit  the 
incorporation  of  certain  signal  generators  which  also  have  application,  such  as  tracking  a 
ramp  or  a  limit  cycle. 

For  lumped  linear  systems,  the  task  of  developing  a  robust  regulator  theory  involves  more 
challenging  cases  when  either: 

a  signal  generator,  or  reference  model,  is  known  but  certain  parameters  in  the  system  or  in 
the  disturbance  channel  are  unknown  {structured  uncertainty)]  or 
a  signal  generator  for  disturbances  is  not  known  (  unstructured  uncertainty) 

Both  of  these  problems  are  central  focal  points  in  the  modern  approach  to  robust  control. 
Setting  some  notation  we  consider  the  problem  of  output  regulation  for  nonlinear  systems 
modeled  by  equations  of  the  form 


x  =  f{x,u,w) 
e  =  h{x,  w), 

with  state  x  e  X  C  R”,  control  input  u  G  R”^,  regulated  output  e  G  R’”  and  exogenous 
disturbance  input  w  eW  CW  generated  by  an  exosystem 

w  =  s{w).  (1-2) 


We  assume  that  f{x,u,w),  h{x,w)  and  w{w)  are  C*  functions  (for  some  large  k)  of  their 
arguments  and  also  that  /(0, 0, 0)  =  0,  h{0, 0)  =  0  and  s(0)  =  0. 

Generally  speaking,  a  problem  of  local  output  regulation  is  to  design  a  feedback  controller 
so  as  to  obtain  a  closed  loop  system  in  which,  when  w{t)  =  0,  a  certain  equilibrium  is  locally 
asymptotically  stable  and,  when  w{t)  ^  0  and  sufficiently  small,  the  regulated  output  e{t) 
asymptotically  decays  to  0  as  t  ^  oo.  The  structure  of  the  controller  usually  depends  on 
the  amount  of  information  available  for  feedback.  We  have  primarily  focused  our  attention 
on  the  case  in  which  the  information  in  question  only  consists,  at  each  time  t,  in  the  value 
e(t)  of  the  error  at  this  time.  In  other  words,  we  consider  controllers  modeled  by  equations 
of  the  form 


^  =  V(^,  e) 

u  =  6(0 


(1.3) 


with  state  ^  G  S°  C  R"",  in  which  77(0  e)  and  6(0  are  functions  of  their  arguments,  and 
7y(0, 0)  =  0,  6(0)  =  0.  The  purpose  of  output  regulation  is  to  obtain  a  closed  loop  system  in 
which,  from  every  initial  condition  in  a  neighborhood  of  the  equilibrium  {x,  0  w)  =  (0, 0, 0), 
the  response  of  the  regulated  output  asymptotically  converges  to  0  as  time  tends  to  00. 
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It  is  important  to  observe  that  output  regulation  entails  the  problems  of  the  asymp¬ 
totic  tracking  of  a  class  of  reference  trajectories,  the  disturbance  attenuation  of  a  class  of 
disturbances,  and  the  requirement  that  both  attenuation  and  tracking  be  achieved  while 
maintaining  internal  stability  of  the  closed-loop  system.  In  this  regard,  one  may  think  of  the 
exosystem  as  consisting  of  two  subsystems,  one  which  generates  signals  to  be  tracked,  and 
one  which  generates  the  disturbances  to  be  attenuated. 

The  case  when  the  exosystem  is  a  harmonic  oscillator  is,  of  course,  classical.  Even  in  this 
special  case,  the  difference  between  state  and  error  measurement  feedback  in  the  problem 
of  output  regulation  is  profound.  To  know  the  initial  condition  of  the  exosystem  is  to  know 
the  amplitude  and  phase  of  the  corresponding  sinusoid.  On  the  other  hand,  to  solve  the 
output  regulation  problem  in  this  case  with  only  error  measurement  feedback  is  to  track,  or 
attenuate,  a  sinusoid  of  known  frequency  but  with  unknown  amplitude  and  phase.  This  is  in 
sharp  contrast  with  alternative  approaches,  such  as  exact  output  tracking,  where  in  lieu  of 
the  assumption  that  a  signal  is  within  a  class  of  signals  generated  by  an  exogenous  system, 
one  instead  assumes  complete  knowledge  of  the  past,  present  and  future  time  history  of  the 
trajectory  to  be  tracked. 

In  this  setup,  the  problem  in  question  can  be  formally  posed  in  the  following  terms. 

Local  Output  Regulation.  Given  a  nonlinear  system  of  the  form  (1.1)  with  exosystem 
(1.2)  find,  if  possible,  a  controller  of  the  form  (1.3)  such  that: 

(a)  the  equilibrium  (x,  0  =  (0)  0)  of  the  unforced  closed  loop  system 

X  =  /(x,  0(0,0) 

^  =  0(^,M^.O)) 

is  locally  asymptotically  stable  in  the  first  approximation, 

(b)  the  forced  closed  loop  system 

X  =  f{x,e{0,w) 

w  =  s{w) 


(1.4) 


(1.5) 


is  such  that 


lim  e(t)  =  0 

t-,00 

for  each  initial  condition  (x(0),O0),  1(^(0))  in  a  neighborhood  of  the  equilibrium  (0,0,0).  < 

The  importance  of  asymptotic  tracking,  disturbance  attenuation,  and  internal  stability 
in  their  own  right  underscores  the  central  role  which  the  problem  of  output  regulation  has 
played  in  the  development  of  classical  and  modern  automatic  control. 

Consider,  for  example,  the  problem  of  asymptotic  output  regulation  for  which  we  observe 
that  we  must  also  be  able  to: 

i.  design  locally  exponentially  stabilizing  state  feedback  laws  for  nonlinear  control  sys¬ 
tems. 
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ii.  determine  conditions  for  the  existence  of  (stable)  forced  oscillations,  in  the  special  case 
when  the  exosystem  is  a  harmonic  oscillator, 

iii.  design  state  feedback  laws  which  shape  the  steady  state  response  to  harmonic  forcing 
for  general  classes  of  exosystems,  as  well  as  for  harmonic  oscillators, 

iv.  design  dynamic  filters,  or  compensators,  which  produce  a  proxy  for  the  system-exosystem 
state  for  feedback  laws  achieving  the  prior  objectives. 

These  are  the  basic  ingredients  to  both  the  output  regulation  problem.  The  key  features 
of  these  basic  ingredients  are  outlined  as  follows: 

Exponential  Stabilization.  This  first  basic  ingredient  is  indeed  a  problem  concerning  the 
linear  approximation 


X  =  Ax  +  Bu 


of  the  nonhnear  control  system 


X  =  f{x,u,  0) 

at  the  equilibrium  {x,u,w)  =  (0,0,0).  In  this  case,  if  all  systems  parameters  are  known, 
problem  (i)  can  be  solved  by  any  of  a  number  of  standard  approaches,  such  as  infinite 
horizon,  linear-quadratic  optimal  control  or  infinite  horizon  Hoo  control  design.  Moreover, 
our  calculations  have  shown  that  the  general  solution  with  error  measurement  feedback 
very  easily  incorporates  any  given  solution  to  problem  (i)  into  a  general  law  achieving  local 
output  regulation.  The  same  calculations  extend  to  the  case  where  certain  plant  parameters, 
denoted  by  /i,  are  unknown.  In  this  case,  one  considers  the  linear  approximation 

X  =  AfiX  +  B^u  (1.6) 


of  a  nonlinear  control  system 

X  =  /^(a:,u,0)  . 

The  problem  of  determining  the  stability  for  arbitrary  families  of  n-dimensional  systems 
with  m-dimensional  inputs  has  been  shown  to  be  NP-hard.  The  problem  of  stabilizing  the 
linearization  (1.6),  for  all  parametric  uncertainties  within  the  class  of  n-dimensional  systems 
with  m-dimensional  inputs,  with  a  linear  controller  is  sometimes  suspected  of  being  to  be 
NP-hard  as  well.  On  the  other  hand,  if  one  admits  nonlinear  controllers  and  insists  that 
the  pair  (A^,  Bfj)  be  controllable  for  every  p,  then  adaptive  stabilization  schemes  are  known, 
underscoring  the  potential  usefulness  of  nonlinear  output  regulation,  even  in  a  linear  context. 

Forced  Oscillations  and  the  Existence  of  a  Steady-State  Response.  For  constant  coeffi¬ 
cient,  linear  control  systems 


X  =  Ax Dw,  w{t)  =  U  sin(a;t) 

where  no  eigenvalue  of  A  lies  on  the  imaginary  axis,  problem  (ii)  has  a  complete  and  satisfying 
resolution.  Viewing  co  as  fixed  by  the  choice  of  the  harmonic  oscillator  with  frequency  uj  as 
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the  exosystem,  for  each  amplitude  U  there  is  exactly  one  initial  condition  x°  which  generates 
a  periodic  trajectory  with  period  T  =  2'k/u.  Moreover,  the  periodic  orbit  is  asymptotically 
stable  if,  and  only  if,  the  unforced  system  is,  i.e.  whenever  the  eigenvalues  of  A  all  lie  in  the 
open  left  half  plane,  which  we  can  assume  has  been  arranged  as  in  problem  (i). 

For  nonlinear  control  systems 


x  =  f(x,0,w)  (1.7) 

where,  for  example,  w{t)  =  U  sin{ut),  the  situation  is  far  more  complex,  with  the  possibility 
of  one,  or  several,  forced  oscillations  with  varying  stability  characteristics  occurring.  In 
addition,  the  fundamental  harmonic  of  these  periodic  responses  may  agree  with  the  frequency 
of  the  forcing  term  (harmonic  oscillations),  or  with  integer  multiples  or  divisors  of  the  forcing 
frequency  (higher  harmonic,  or  subhaxmonic,  oscillations).  Despite  a  vast  literature  on 
nonlinear  oscillations,  a  subject  with  its  origins  in  celestial  mechanics,  only  for  second  order 
systems  is  there  much  known  about  the  stability  of  forced  oscillation  and,  in  particular, 
which  of  these  three  kinds  of  periodic  responses  might  be  asymptotically  stable. 

It  is  interesting,  then,  to  view  the  solution  of  the  problem  of  output  regulation  for  local 
and  more  global  problems  in  this  context.  We  shall  limit  ourselves  to  a  discussion  about  the 
local  case.  Suppose,  for  simplicity  of  discussion,  that  (1.7)  is  affine  in  w,  i.e.,  so  that  we  have 

^  =  fo{x)+Pix)w  , 


or,  in  particular. 


X  =  fo{x)  +p{x)U  sin(a;t) 

and  we  can  view  p{x)U sin(a;t)  as  a  small  perturbation.  This  is  the  basis  of  the  method  of 
averaging  for  determining  the  existence  and  stability  of  periodic  orbits.  A  similar  analysis 
is  motivation  for  the  method  of  “harmonic  balance”  in  which  a  Fourier  series  is  considered 
for  an  assumed  periodic  trajectory  of  period  T  =  2'k/u  and  the  Fourier  coefficients  are 
determined  so  that  the  differential  equation  is  satisfied;  i.e.,  so  that 

f{x{s),0,  U 8m{ojs))ds  =  0  .  (1.8) 

Averaging  and  harmonic  balance  have  their  origin  in  the  method  of  “small  parameters,” 
pioneered  by  Poincare.  Briefly,  in  the  case  of  constant  coefficients  linear  control  systems, 
the  initial  condition  x°  which  generates  a  periodic  trajectory  can  be  expressed  as  a  linear 
function  of  the  amplitude  U 


x°  =  PU 


using  the  variation  of  parameters  formula,  so  that 


P  =  {I- 


qMx  ®)_Dsin(a;s)(is  . 
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Poincare’s  idea  was  to  take  the  linear  function  obtained  for  the  first  approximation  as  the 
first  term  in  a  power  series  representation  in  U, 

x°  =  PU  +  P2C/2  +  PsU^  +  •  •  • 

which  can  in  principle  be  developed  from  (1.8)  for  an  assumed  periodic  trajectory  of  period 
T  or,  alternatively,  for  an  initial  condition  which  generates  this  trajectory.  Of  course,  one 
needs  to  check  that  the  series  converges  and  represents  a  function  7c{U)  in  a  neighborhood 
of  C/  =  0.  Our  approach  is  geometric  and  seeks  instead  to  “sum  the  series” 

=  PU  +  PiU^  +  PsU^  +  ...  =  Tr{U) 

by  characterizing  tt{U)  as  the  solution  of  a  system  of  partial  differential  equations,  for  which 
the  center  manifold  theorem  guarantees  existence  of  a  solution.  More  explicitly,  observe  that 
in  general  the  response  of  system  (1-7),  in  the  initial  state  a;(0)  =  xo  and  subject  to  the  input 
U  sm{ut),  coincides  with  the  response  of  the  autonomous  system 

x  =  f{x,0,w))  qs 

w^Sw  ^ 


where  w  G  K.^, 

5=  f  ^ 

—u 

and 

x(0)  =  Xo  , 

Assuming  that  we  have  solved,  and  implemented,  the  feedback  stabilization  problem,  the 
center  manifold  can  be  expressed  as  the  graph  of  a  mapping  x  =  ^{w),  where  Tr{w)  is  a 
function  satisfying 

^Sw  =  f{-K{w),d{w)) 


u 


,  d{w)  =  Wi  , 


w(0)  =  1^0  = 


and  7r(0)  =  0. 

The  restriction  of  the  flow  of  (1.9)  to  its  center  manifold  is  indeed  a  copy  of  the  flow  of 


w  =  Sw  . 


As  a  consequence,  if  U  is  sufficiently  small  and  Xq  =  7r(t(;o),  the  response  x{t)  exists  for 
all  t  G  E  and  is  periodic  of  period  T.  Moreover,  once  one  has  implemented  a  solution  to 
problem  (i),  asymptotic  (orbital)  stability  of  the  steady-state  response  also  follows  from  the 
center  manifold  theorem.  Furthermore,  for  the  problem  of  output  regulation  treated  in  our 
work,  the  center  manifold  can  be  shown  to  be  unique. 


6 


For  a  nonlinear  control  system  with  unknown  parameters 

X  =  u,  0) 

these  arguments  persist,  yielding  an  asymptotically  stable,  steady-state  response  which,  of 
course,  is  dependent  on  //,  and  hence  unknown  for  purposes  of  problem  (iii). 

Shaping  the  Steady-State  Response.  The  third  basic  ingredient  in  output  regulation,  the 
ability  to  shape  the  steady-state  response  of  a  nonlinear  control  system 

X  =  f{x,u,w) 

is  prescriptive  in  nature,  and  would  appear  to  be  quite  challenging  from  the  point  of  view 
of  the  more  descriptive  tools  for  predicting  the  existence  of  a  steady-state  response.  In  our 
approach  to  classical  and  robust  output  regulation,  the  existence  of  a  steady-state  response 
and  the  ability  to  shape  this  response  are  embodied  in  each  of  two  equations,  known  as  the 
“regulator  equations,” 


dir 

dw 


s{w)  — 


f{n{w),c{w),w) 
h{n{w),  w) 


(1.10) 


for  two  C'^  (with  A:  >  1)  mappings  tt  :  Wo  — >  R”  and  c  :  Wq  — >•  (where  Wo  C  W  is  a 

neighborhood  of  w  =  0),  with  7r(0)  =  0  and  c(0)  =  0,  and  whose  solution  necessarily  exists 
if  the  problem  of  local  output  regulation  can  be  solved  using  any  method. 

While  the  first  of  the  regulator  equations  will  always  have  a  solution,  as  a  consequence 
of  the  center  manifold  theorem,  the  solvability  of  the  system  of  regulator  equations  can  be 
expressed  in  a  system  theoretic  framework  which  has  an  appealing  “frequency  domain”  inter¬ 
pretation.  In  classical  automatic  control,  there  is  a  simple,  intuitive  condition  for  solvability 
of  the  problem  of  output  regulation:  no  transmission  zero  of  the  plant  to  be  controlled  should 
coincide  with  a  natmal  frequency  of  the  signal  to  be  tracked  (or  attenuated).  For,  while  a 
unique  steady-state  periodic  response  to  harmonic  forcing  at  such  a  frequency  certainly  exists 
for  a  linear,  single-input  single-output  transfer  function,  if  the  emitted  response  is  absorbed 
at  this  frequency  we  cannot  adjust  its  amplitude  or  phase  by  feedback. 

This  fundamental  condition  can  be  recast  in  a  state-space  form  by  introducing  a  linear 
operator  whose  spectrum  on  one  subspace  coincides  with  the  plant  transmission  zeros  and  on 
whose  spectrum  on  another  subspace  coincides  with  the  natural  frequencies  of  the  exosystem. 
In  this  setting,  to  say  that  no  transmission  zero  of  the  plant  coincides  with  a  natural  frequency 
of  the  signals  to  be  tracked  is  to  say  that  these  two  subspaces  are  complementary.  We  can 
also  describe  an  existence  theory  for  the  regulator  equations  which  provides  a  nonlinear 
enhancement  of  this  criterion,  in  terms  of  zero  dynamics.  Briefly,  since  whenever  w  =  0 
and  the  system  output  h{x,  0)  is  constrained  to  be  zero  we  must  have  that  the  error  is  zero, 
the  zero  dynamics  of  the  augmented  system  contains  the  zero  dynamics  of  the  system  to  be 
controlled.  In  this  setting,  the  regulator  equations  are  solvable  just  in  case  the  plant  zero 
dynamics  are  complemented  in  the  augmented  zero  dynamics  by  a  copy  of  the  exosystem. 
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Filtering  a  Proxy  for  the  Plant/ Exosystem  State.  If  the  plant/exo-system  is  exponen¬ 
tially  detectable  using  the  error  variables  as  an  output  measurement,  one  can  verify  that  a 
“separation  principle”  holds  for  output  regulation  via  state  feedback  and  a  standard  state 
observer  scheme.^  Unlike  exponential  stabilizability  of  the  plant,  however,  detectability  of 
the  augmented  system  is  not  a  necessary  condition  for  the  solvability  of  the  problem  of  out¬ 
put  regulation.  More  importantly,  detectability  of  the  augmented  system  will  not  hold  when 
we  treat  the  case  of  robustness  with  respect  to  real,  parametric  uncertainty.  Flather,  the 
derivation  of  the  regulator  equations  for  the  error  measurement  case  reveals  that  the  actual 
necessary  condition  relating  to  measurements  of  the  augmented  state  is,  in  fact,  a  geometric 
formulation  of  the  “internal  model  principle.”  In  short,  rather  than  it  being  necessary  to 
be  able  to  recover  the  augmented  state,  it  is  only  necessary  that  any  dynamic  compensator 
which  achieves  output  regulation  must  also  contain  a  copy  of  the  exosystem. 

With  these  motivations  in  mind,  we  now  combine  the  actual  necessary  conditions,  em¬ 
bodied  in  the  regulator  equations,  with  the  notion  of  “immersion”  of  a  nonlinear  system  to 
derive  a  new  class  of  systems  which  generalize  dynamic  observers  and  provide  a  “proxy”  for 
the  state  in  a  compensator  design  which  provides  an  alternative  to  the  design  based  on  a 
separation  principle. 

Consider  the  problem  of  designing  feedback  laws  achieving  output  regulation  for  periodic 
signals  with  an  unknown  frequency.  For  simplicity,  we  first  consider  a  one-dimensional  scalar 
linear  control  system. 

For  a  linear  plant: 


X  =  u 

y  =  h{x)  =  X 


and  an  exogenous  system  defined  by  the  harmonic  oscillator  w  =  s{w) 


SWi  =  W2 
W2  = 
fi  =  0 
■k{w)  =  Wi 

with  an  unknown  by  fixed  frequency.  We  first  note  that  the  harmonic  oscillator  is  now  rep¬ 
resented  by  a  (critical)  nonlinear  system.  We  wish  to  design  a  (possibly  nonlinear)  feedback 
system  so  that  the  plant  state  x  tracks  tt{w)  for  a  continuous  range  fx  G  [.5, 4],  In  this  case, 
we  define  the  tracking  error  to  be  e  =  a;  —  'k{w)  and  proceed  to  formulate  and  solve  the 
regulator  equations,  which  of  course  are  quite  simple  in  this  case. 

The  second  regulator  equation  implies 


X  =  'k{w)  =  Wi 


which,  when  differentiated  yields  the  first  regulator  equation  and  an  explicit  expression  for 
the  zero  error  input,  u  —  c{w). 

^See  [12]. 
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Now  define  a  feedback  law 


dTT{w)  . 

X  =  - - S{w)  =  W2, 

ow 

c{w)  =  W2{  zero  error  input 


} 


u  =  c{w)  +  k{x  —  wi) 

where  k  is  chosen  to  stabilize  the  plant  in  the  first  approximation  when  7r(w)  =  0.  The 
closed  loop  system  is  therefore  given  by 

X  =  W2  +  k{x  —  Wi) 

Wl  =  W2 
W2  =  -/x^wi 
fi  =  0 
e  —  x  —  wi 

leading  to  the  following  performance. 


mu  0  .  mu  0 


Figure  2.2.1  Figure  2.2.2 

We  note,  however,  that  although  the  feedback  scheme  does  not  explicitly  contain  the 
unknown  parameter  fi  it  does  contain  the  exosystem  state  variable  W2  which  is  not  observable 
from  the  error  without  explicit  knowledge  of  /i.  Actually,  we  need  only  produce  the  desired 
feedback  law  asymptotically,  not  exactly  at  each  instant  of  time,  in  order  to  solve  this  robust 
output  regulation  problem.  That  is,  in  lieu  of  an  observer  based  design,  we  need  only  design 
a  dynamical  system  which  asymptotically  produces  the  error  zeroing  control. 

Next  we  illustrate  this  design  philosophy  for  a  periodic  motion  with  unknown  frequency 
and  phase  but  with  known  amplitude;  i.e.,  we  assume  that  the  norm  of  the  initial  condition 
w(0)  is  known  and,  for  simplicity  is  equal  to  1. 

In  fact,  robust  tracking  problem  is  intimately  related  to  a  problem  of  output  regulation 
for  a  four-dimensional  nonlinear  oscillator  which  produces  oscillations  at  every  frequency.  In 
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fact,  this  system  was  designed  to  asymptotically  produce  the  control  law  designed  above. 
More  explicitly,  consider  the  harmonic  oscillator  in  its  more  typical  representation 


=  fJ,W2 
W2  =  -fJ>Wi 

fi,  =  0 

'k{w)  =  Wi 

Suppose  one  now  defines  =  Wi  and  begins  differentiating  m  times  until  •••Cm-i)- 

One  then  has  a  closed  system  whose  linear  approximation  is  observable  and  whose  first  state 
=  Wi  for  the  “right”  initial  conditions. 

il=Wi=  HW2  =  6 
^2  =  -j/wi  =  ^3 
^3  =  —fJ'^W2  =  ^4 
^4  =  =  ^5 

If  the  restriction  taf  +  tai  =  1  is  imposed,  then 

^4  =  +  wl) 


which  can  be  written  as 

4  =  6^1  - 

Now  if  wi(0)  =  Ui  and  1^2(0)  =  02  then  the  four-dimensional  nonlinear  system  simulates 
the  two-dimensional  system  for  the  “right”  conditions  stated  below. 


6  =  ^2  1 

'  6(0)  =  ai  ^ 

^2  =  ^3 

6(0)  =  na2 

6  =  ^4 

>  subject  to  < 

6(0)  =  > 

^4  =  6 

6(0)  =  -)a^02 

aj  +  al  =  l  J 

Furthermore,  this  same  four-dimensional  nonlinear  system  oscillates  at  a  frequency  which 
is  determined  by  its  initial  conditions.  Finally,  we  note  that  this  system,  with  the  first  coor¬ 
dinate  as  the  system  output,  is  in  fact  observable.  In  particular,  one  can  recover  functionals 
of  the  state,  such  as  th  second  coordinate,  and  so  the  nonlinear  oscillator  can  be  used  as  a 
proxy  for  the  feedback  law  designed  for  the  harmonic  oscillator  having  arbitrary  frequency. 

Briefly,  the  second  regulator  equation  yields  7r(^)  =  which,  when  differentiated,  yields 
the  first  regulator  equation  and  an  explicit  expression  for  the  zero  error  input,  u  =  c{^) 


c(0  =  '^2  {zero  error  input}. 

Now  define  a  feedback  law 

u  =  c{^)  +  k{x  -  ^i) 

where  the  scalar  k  is  chosen  to  stabilize  the  one-dimensional  plant.  The  closed  loop  system 
is  therefore  given  by 

x  =  ^2  +  k{^-  ^i) 

=  6 
6  =  6 
6  =  ^4 

^4  =  6^3  ”  ^26 
e  =  x-^i 


Notice  that  in  reality  the  task  has  not  yet  been  truly  completed  because  the  state  ^2 
is  not  actually  available  for  measurement.  However,  for  a  neutrally  stable  system  having 
a  detectable  linear  approximation,  there  exists  a  local  exponential  observer  [1].  Recall  the 
nonlinear  exosystem  with  output  7(^) 


=  6 
6  =  6 
6  =  ^4 

^4  =  ”  ^26 

7(^)  =  6 


The  linear  approximation  is  shown  below 


A  = 


0  10  0 
0  0  10 
0  0  0  1’ 
0  0  0  0 


(7=  [  1  0  0  0  ]  . 


It  is  not  difficult  to  show  that 

C 

CA 

CA^ 

CA^ 


has  full  rank.  In  accordance  with  [1],  let  eo  =  2:1  —  be  the  error  between  the  first  observer 
state  and  the  only  available  state  from  the  exosystem  and  seek  L  such  that  A  —  LC  is 
Hurwitz.  One  such  choice  is 


10 

35 

50 

24 
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The  resulting  nonlinear  observer  is 


— 1 

o 

o 

o 

_ 1 

■  10  ■ 

0 

0  0  10 

35 

0 

0  0  0  1 

+ 

50 

6o  + 

0 

0  0  0  0 

24 

.  Zizj  -  Z^Zs  _ 

Thus,  the  state  estimates  the  inaccessible  state  ^2  in  the  feedback  loop. 

We  now  depict  the  plots  of  the  tracking  state  Xi  and  the  tracking  error  e  =  x-^i  along 
with  the  observed  feedback  state  Z2  and  the  observing  error  62  =  Z2  —  6  nre  shown. 


plot  of  x1(t)  for  mu  =  .5  to  2 


mu  0.5  0  tfrom0to3*pi 
plot  of  z2(t)  for  mu  =  .5  to  2 


plot  of  e1  (t)  for  mu  =  .5  to  2 


0  tfrom0to3*pl 
plot  of  e2(t)  for  mu  =  .5  to  2 


mu  0.5  0  t  from  0  to  3*pi  ^  t  from  0  to  3*pi 


Figure  2.2.3 


We  shall  now  describe,  in  a  more  formal  and  complete  manner,  the  concepts,  construc¬ 
tions,  and  basic  theory  underlying  the  use  of  “asymptotic  proxies”  in  feedback  design  for 
unobservable,  or  nondetectable,  systems  and  its  role  in  robust  control.  As  it  turns  out,  a  key 
construction  to  understanding  this  technique  is  the  formulation  of  the  problem  of  system 
immersion,  pioneered  by  M.  Fliess.  As  already  we  have  already  observed,  the  first  one  of  the 
regulator  equations  (1.10)  expresses  the  property  that  the  graph  of  the  mapping  x  =  'ir{w) 
is  an  invariant  manifold  for  the  composite  system 

x  =  f{x,c{w),w) 
w  =  s{w) 


while  the  second  equation  expresses  the  property  that  the  error  map  e  =  h{x,  w)  is  zero  at 
each  point  of  this  invariant  manifold.  In  slightly  different  (and  indeed  equivalent)  terms, 
one  may  say  that  the  regulator  equations  call  for  the  existence  of  a  map  u  =  c{w)  with 


12 
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the  property  that  the  composite  system  (1-11),  namely  the  controlled  plant  driven  by  the 
autonomous  subsystem 


w  =  s{w) 
u  —  c{w) 


(1.12) 


has  an  invariant  manifold  on  which  the  error  map  is  zero. 

The  existence  of  such  an  invariant  manifold,  however,  is  not  the  only  condition  that 
system  (1.11)  is  required  to  satisfy.  In  fact,  by  extending  the  analysis  presented  in  [12], 
it  is  possible  to  identify  additional  necessary  conditions  (which  are  very  similar  to  those 
established  earlier  in  [1],  [2])  for  the  solvability  of  a  problem  of  output  regulation.  These 
additional  conditions  can  be  expressed  as  properties  of  the  autonomous  subsystem  (1.12)  and, 
together  with  those  expressed  in  (1.10),  can  be  used  to  provide  a  complete  set  of  necessary 
and  sufficient  for  the  solvability  of  the  problem  under  investigation. 

The  formulation  of  these  new  conditions  is  based  on  the  notion  of  immersion  of  a  given 
system  into  another  one,  which  is  defined  as  follows.  Let  {X,  /,  h}  denote  the  autonomous 
system 


X  =  f{x) 

y  =  h(a:). 


(1.13) 


with  state  x  E  X  and  output  y  G  in  which  we  suppose  /  to  be  a  C'‘  vector  field  and  h 
a  mapping,  for  some  sufficiently  large  k,  with  /(O)  =  0  and  h(0)  =  0. 

Definition  1.1.  (see  [3],  [4])-  System  {X,  f,  h}  is  immersed  into  system  {X',  f,  h'}  if  there 
exists  a  C’‘  a  mapping  r  :  X  — >  X',  satisfying  r(0)  =  0  and 


h{x)  h{z)  ^  h'{r{x))  ^  h'(r(^)). 


such  that 

|^/(a^)  =  f{r{x)) 
h{x)  —  /i'(r(x)) 

for  all  X  E  X. 


Note  that,  if  {X,  /,  h}  is  immersed  into  {X',  /',  h'},  any  output  response  of  {X,  /,  h}  (i.e. 
the  response  from  any  arbitrary  initial  state  xq  E  X)  is  also  an  output  response  of  {X',  f,  h'} 
(namely,  the  response  from  the  initial  state  r(xo)  E  X'). 

For  example  the  system 


Wl  =  pLW2 
W2  =  -p-Wi 
/{  =  0 
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evolving  on  the  state  manifolds 

S^{wi,W2,lj)  :  +W2  =  1} 

can  be  immersed  into  the  4-dimensional  nonlinear  oscillator  via  the  mapping 


■^r 

r  *1 

m 

C2 

Wi 

IJLW2 

6 

=  T 

W2 

-iBwi 

-[x'^W2_ 

Thus,  the  nonlinear  oscillator  contains  an  “internal  model”  for  unit  amplitude  sinusoidal 
motion  of  arbitrary  phase  and  period,  a  property  which  enables  it  to  asymptotically  repro¬ 
duce  the  kind  of  feedback  laws  needed  to  shape  the  steady-state  response  of  the  closed-loop 
system  designed  above.  The  methods  used  in  this  example  actually  extend  quite  far,  under¬ 
scoring  the  utility  of  these  concepts  for  the  solvability  of  output  regulation  problems  with 
real  parametric  uncertainty. 

(i)  there  exist  (7^  (with  A:  >  1)  solutions  tt  :  Wo  — >  3?”  and  c  :  Wo  9?"^  (where  Wq  C  W  is 
a  neighborhood  of  w  =  0),  with  7r(0)  =  0  and  c(0)  =  0,  of  the  regulator  equations  (1.10) 

(ii)  there  exist  an  integer  u,  a  neighborhood  Ho  of  the  origin  in  S?*',  a  vector  field  (p  : 
So  — >  TEq  and  a  mapping  7  :  So  — >  5?”^,  with  (/?(0)  =  0  and  7(0)  =  0,  such  that 


i  =  ^(0 

u  =  7(0 


(1.14) 


is  locally  exponentially  detectable  at  ^  =  0  and  {Wo,5,  c}  is  immersed  into  {So,</?,7}) 

(iii)  the  linear  approximation  of 

X=|(0,0,0).  B  =  |(0.0,0).  c=|(0.0). 

at  the  equilibrium  {x,  u,  w)  =  (0, 0, 0)  is  stabilizable  and  detectable,  and  the  matrix 


A -XI  B 
C  0 


is  nonsingular  for  each  eigenvalue  A  of 

o  =  |(0) 

having  nonnegative  real  part. 

As  a  special  case,  consider  a  system  modeled  by  equations  of  the  form 

X  =  f{x,u,w,fj) 
e  =  h{x,w,fj,) 


(1.15) 


(1.16) 
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with  state  x  6  R”,  control  input  ueW^,  regulated  output  e  G  subject  to  an  exogenous 
disturbance  input  w  G  R'^,  in  which  /x  G  T’  C  R^  is  a  vector  of  unknown  parameters  and  V  is 
compact  set.  f{x,  u,  w,  ix)  and  h{x,  w,  fi)  are  functions  of  their  arguments  (for  some  large 
k),  and  /(0, 0, 0,  /r)  =  0  and  h{0, 0,fj,)=0  for  each  value  of  fx.  Without  loss  of  generality,  we 
suppose  0  G  int(P).  The  exosystem 


w  =  Sw 

(1.17) 

is  assumed  to  be  neutrally  stable. 

The  control  input  to  (1.16)  is  to  be  provided  by  a 
form 

controller  modeled  by  equations  of  the 

i  = 

u  =  e{() 

(1.18) 

with  state  ^  G  R'',  in  which  ri{^,e)  and  0{^)  are  C’^  functions  of  their  arguments,  and 
77(0,0)  =  0,  0(0)  =  0. 

In  order  to  motivate  some  of  our  subsequent  hypotheses,  we  recall  first  a  necessary  condi¬ 
tion  for  the  solution  of  the  problem  of  structurally  stable  output  regulation,  which  is  indeed 
a  prerequisite  for  robust  regulation. 

Structurally  stable  local  output  regulation.  Given  a  nonlinear  system  of  the  form  (1.16) 
with  exosystem 

w  =  s{w)  , 

(1.19) 

find  a  controller  of  the  form 

II  II 

(1.20) 

such  that,  for  some  neighborhood  "P  of  /i  =  0  in  R^  and  for  each  ^  G  P: 

(a)  the  equilibrium  (x,  =  (0, 0)  of  the  unforced  closed  loop  system 

X  =  f{x,e{^),0,H) 

C  =  ix)) 

is  locally  asymptotically  stable  in  the  first  approximation, 

(b)  the  forced  closed  loop  system 

X  =  f{x,e{^),w,fx) 
i  =  'n{i,h{x,w,ix)) 
w  =  s{w) 

is  such  that 

lim  e{t)  =  0 


15 


for  each  initial  condition  (x(0),^(0),tt;(0))  in  a  neighborhood  of  the  equilibrium  (0,0,0). 

The  following  result  provides  necessary  and  sufl&cient  conditions  for  the  existence  of 
solutions  of  the  problem  structurally  stable  output  regulation.  Set 


A{fj)  = 


df 

dx 


B(ii)  = 


J  (0,0, 0,m) 


du 


J  (0,0,0,;i) 


dh 

dx 


(0,0,/i) 


Theorem  1.1.  Suppose  the  exosystem  (1-17)  is  neutrally  stable.  There  exists  a  solution 
of  the  problem  of  structurally  stable  output  regulation  if  and  only  if  there  exist  mappings 
X  =  ir^{w,pL)  and  u  =  c^{w,p,),  with  7r®^(0, /r)  =  0  and  0^(0,//)  =  0,  both  defined  in  a 
neighborhood  W°  xV°cMf'xVof  the  origin,  satisfying  the  conditions 

dw'  ""  f{T^^{w,(J-),c^{w,fa),w,iJ,)  ^^  21) 

0  =  h('K^{w,  p,),w,  n) 

for  all  {w,  p)  €  W°  X  V° ,  and  such  that  the  autonomous  system  with  output 

w  =  Sw 

p  =  0  (1.22) 

u  =  c^{w,  p) 


is  immersed  into  a  system 


i  =  <P{0 
u  =  7(0  . 


defined  on  a  neighborhood  S°  of  the  origin  in  R'',  in  which  (/:(0)  =  0  and  7(0)  =  0  and  the 
two  matrices 


d(p 

,  r  = 

57 

(=0 

are  such  that  the  pair 


■  ^4(0)  0' 

'B(0) 

_NC{0)  $ 

J 

0 

(1.23) 


is  stabilizable  for  some  choice  of  the  matrix  N,  and  the  pair 


[(7(0)  0], 


■a(o)  jB(o)r 

0  $ 


(1.24) 


is  detectable. 


Remark  1.1.  Note  that  the  linear  approximation  of  (1.22)  at  the  equilibrium  (w,  p)  =  (0, 0) 
cannot  be  detectable.  In  fact,  since  c®'(0,  p)  =0  by  hypothesis. 


dp 


(0,  p)  =  0  , 
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and  the  linear  approximation  in  question  is  characterized  by  a  pair  of  matrices  of  the  form 

S  O' 

0  0 

which  is  indeed  not  detectable.  Thus,  it  is  not  possible  to  have  the  conditions  of  the  Theorem 
directly  satisfied  by  the  trivial  immersion  of  (1.22)  into  itself.  However,  as  shown  below, 
(1.22)  may  be  immersed  into  another  system,  having  a  detectable  approximation  at  ^  =  0. 

Remark  1.2.  Note  that  (1.22)  is  immersed  into  a  (possibly  nonlinear)  system  having  an 
observable  linear  approximation  if  there  exists  a  function  4>{Co,  Ci)  •  •  •  >  C9-1)  that 

Llc^iw,  p)  =  (l){c^{w,  p),  LgC^iw,  p),...  ,  Ll~^c^{w,  p))  . 

In  particular  (see  [3]),  system  (1.22)  is  immersed  into  a  linear  system  if  and  only  if,  for 
some  set  of  real  numbers  oq,  Oi, . . .  ,  Og-i, 

V^c^iw,  p)  =  aoc^{w,  p)  +  aiLsC^{w,  /r)  +  •  •  •  +  ag-iLg  ^c^{w,  p)  ,  (1-25) 

for  all  {w,p)  E  W°  X  V°. 

1.2  Boundary  Control  of  Distributed  Parameter  Systems 

One  of  our  longer  term  goals  studied  in  this  research  project  was  the  development  of  a 
systematic  feedback  design  methodology  for  distributed  parameter  systems  which  would 
retain  some  of  the  intuitive  appeal  of  classical  automatic  control  in  a  spirit  similar  to  the 
program  pursued  for  the  control  of  lumped  systems. 

This  program  included  the  development  of  the  concept  of  zero  dynamics  for  nonlinear  dis¬ 
tributed  parameter  systems  and  of  a  feedback  design  theory  based  on  asymptotic  properties 
of  the  zero  dynamics.  We  have  approached  this  problem  based  in  part  on  the  development 
of  a  root-locus  methodology  for  distributed  parameter  systems  given  in  [93] . 

In  particular,  for  a  very  general  class  of  parabolic  boundary  control  systems  we  were  able 
to  present  a  fairly  complete  picture  of  root-locus  design  methodology  based  on  analysis  of 
the  closed  loop  transfer  function.  Namely,  we  were  able  to  identify  the  open  loop  zeros  and 
poles  of  the  meromorphic  open  loop  transfer  function.  Then  by  introducing  a  proportional 
error  boundary  feedback  control  law  the  resulting  closed  loop  poles  were  determined  from 
the  return  difference  equation.  For  a  single  input  single  output  system,  we  were  able  to 
show  that  the  closed  loop  poles  varied  from  the  open  loop  poles  to  the  open  loop  zeros 
as  the  gain  parameter  varied  from  zero  to  plus  or  minus  infinity  depending  on  the  sign  of 
the  instantaneous  gain.  This  analysis  provided  a  means  for  establishing  boundary  feedback 
stabilization  for  minimum  phase  systems  based  on  fairly  simple  root-locus  design  methods. 

Based  on  this  work  we  were  able  to  identify  a  preliminary  definition  of  zero  dynamics 
for  boundary  control  systems  which  provides  an  explicit  representation  of  the  zero  dynamics 
as  a  distributed  system  whose  dynamics  are  governed  by  the  zeros  of  the  open  loop  system. 
Setting  notation,  consider  a  general  linear  or  nonlinear  distributed  parameter  system 

w  =  A{w)  (1.26) 


[»  0]. 
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with  inputs  and  outputs 

u(t)  =  B{w){t),  y{t)  =  C{w){t) 

and  initial  conditions 

ii;(0)  =  Wo- 

Associated  with  this  system  we  define  the  zero  dynamics  to  be  the  system  with  the  out¬ 
put  constrained  to  be  zero,  i.e.,  C{w){t)  =  0  In  particular,  if  Wo  satisfies  the  constraint 
C{wo){0)  =  0  and  w{t)  is  a  solution  of  the  zero  dynamics,  then  u{t)  =  B{w){t)  is  a  control 
which  “maintains”  the  constraint  y{t)  =  0. 

This  can  be  viewed  within  the  context  of  root-locus  as  follows.  We  introduce  a  propor¬ 
tional  error  feedback  in  the  form 


u  =  —ky,  k  eR. 


This  provides  a  closed  loop  system 


w  =  A{w)  (1-27) 

B{w){t)  +  kC{w){t)  =  0 

w;(0)  =  Wo- 


For  this  closed  loop  system,  we  see  that  the  uncontrolled  problem  is  obtained  when  k  =  0 
and  in  the  high  gain  limit,  at  least  at  a  formal  level  we  obtain  the  zero  dynamics  as  k  tends 
to  infinity 


w  =  A{w)  (1-28) 

C{w){t)  =  0 

tt;(0)  =  Wo- 

This  turns  out  to  be  a  very  natural  definition  of  zero  dynamics  within  the  context  of 
boundary  feedback  control  problems  for  distributed  parameter  systems  governed  by  partial 
differential  equations.  As  in  the  finite  dimensional  case,  we  have  been  able  to  prove,  under 
appropriate  hypotheses,  that  as  we  tune  a  gain  parameter  k,  the  closed-loop  trajectories 
approach  the  trajectories  of  the  zero  dynamics.  For  linear  systems  we  have  been  able  to 
prove  this  result  based  on  our  root-locus  analysis  in  [93].  For  nonlinear  systems,  in  particular 
for  a  boundary  controlled  Burgers’  equation,  we  have  been  able  to  prove  a  result  of  this 
type  in  [89,  90]  using  semigroup  methods.  These  results  required  that  the  initial  data  be 
sufficently  small  and  smooth.  More  recently,  based  on  our  work  in  [83],  we  have  been  able  to 
extend  these  result  to  include  arbitrary  initial  data  using  energy  methods  and  Galerkin 
approximations.  This  work  is  currently  in  press.  Based  on  this  analysis  we  have  been  able  to 
carryout  significant  numerical  simulations  and  testing  as  proposed  in  our  research  objectives 
for  this  effort. 
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Setting  notation,  the  forced  and  controlled  Burgers’  equation  on  the  finite  interval  (0, 1) 
is 

Wt  -  ewxx  +  {f{'w))x  =  ^)> 

w{x,  0)  =  (p{x),  (p  €  L^{0, 1) 


where  f{vj)  =  — . 

We  introduce  a  feedback  mechanism  by  the  law: 

(/(<)  =  - (*2  «e[0,oo). 

where  ko  >  0,ki  >  0  are  the  gain  parameters.  This  feedback  law  can  be  rewritten  in  the 
form  of  “radiation”  boundary  conditions 

Wx{0,t)  -  kow{Q,t)  =  0 
Wx{l,  t)  +  ki  io(l,  t)  =  0. 

The  zero  dynamics  is  defined  by  the  condition  that  the  output  be  constrained  to  be  zero: 

y(t)  =  0  ,  t  G  [0,  oo). 

This  condition,  which  formally  corresponds  to  passing  ko  and  ki  to  infinity,  gives  the  Dirichlet 
boundary  conditions 

iy(0,  t)  —  w{l,  t)  =  0  ,  t  G  [0,  oo). 

Closed  Loop  System 

Wt  -  ewxx  +  {f{w))x  =  h(x,  t),  (1.29) 

w{x,  0)  =  ip  G  L^(0, 1) 

—  Wx{0,  t)  +  kow{0,  t)  —0 
Wx(l,  t)  +  kiw{l,  t)  =  0 

and 

The  Zero  Dynamics 
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Wt  -  eWxx  +  {f{'w))x  =  h{x,  t), 

w{x,  0)  =  (p{x),  (p  G  1/^(0, 1) 
w{0,t)  =  0 
w{l,t)  —  0 

Let  w'‘{x,  t)  denote  the  solution  of  the  closed  loop  boundary  control  problem  corresponding 
to  {ko,ki)  and  w{x,t)  of  the  zero  dynamics  problem. 

Assume  that  (p  and  h  in  both  problems  are  the  same.  Then 

sup  \\w\t)  -  w{t)\\  0. 

te[o,T] 


If  in  addition,  the  forcing  term  h  satisfies  the  extra  condition: 

hf  G  ([0, 1]  X  [to,  T]) 

then  we  can  show 

sup  ||w^(t)  -  ty(<)|liri(o,i)  0 

t€[fo,r] 

for  all  to  >  0. 

Due  to  the  embedding  H^{0, 1)  C  C[0,  l],this  implies  that  for  t  >  0 

max  sup  |u;^(a;,t)  —  w(x,t)|  °°>  0. 

Commentary  on  the  proof  of  convergence: 

We  have 

+  f{w%  =  h{x) 
w^{0,t)  —  kQw’‘{0,t)  =  0 
w^(l,t)  +  kiw'‘{l,t)  =  0 
w'^{x,0)  =  (f>{x), 

and 

Wt  -  +  fiw\  =  h{x) 

w‘^{0,t)  =  0 
=  0 
0)  = 
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We  let 

u{x,  t)  =  w’^{x,  t)  —  w^{x,  t) 

and  show  that 

lim  \\u{t)\\  =  0,  uniformly  in  i  G  [0,r]. 

It  follows  that  u  satisfies 

Ut  -  +  ifiw'')  -  f{w°))x 

Ux(0,t)  -  kou{0,t)  =  -wl{Q,t) 

+  kou{l,t)  = 


\j^\\nf  +  AW\\l  =  -  jjnw^)  -  f(w°)Udx 

+£f(«), 

where 

F(t)  =  u{0,t)w°{0,t) 
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Now  using  a  variety  of  estimates  (Poincare,  Cauchy-Schwartz,  Young,  Sobolev,  Gronwall, 
etc),  we  estimate  all  the  terms  on  the  right  and  arrive  at  an  expression 


+  C2{to,T)\\u{to)f 

for  all  0  <  to  <  T  <  oo.  Here  Ci{to,T)  is  a  continous  function  of  to  >  0  ^■iid  behaves  like 
to  near  to  =  0.  C2  is  continous  in  to  and  bounded  at  to  =  0. 

Consider  the  following  statements: 

1.  ||M(t)|l  <  |lw*(t)  -  (f)\\  +  ||w°(t)  -  (j)\\. 

2.  w’^{t)  (f)  in  L^(0, 1)  uniformly  with  respect  to  ko  and  ki. 

3.  The  above  convergence  is  uniform  with  respect  to  ko  and  ki  if  ko,ki  G  [k,  00),  where 
k>0. 

4.  Similarly,  ||w°(t)  —  (j)\\  0. 

5.  w*(t)  and  w^{t)  are  continuous  functions  of  t  in  L^-norm. 

6.  It  follows  from  parts  1.-  4.  that 

ii«(t„)ii  ^  0 

uniformly  with  respect  to  ko  and  ki,  if 

ko,  ki  G  [k,  00),  k  >  0. 


Take  any  (5>  0,  then  from  statement  6.  above,  it  follows  that  we  can  select  to  G  (0,  T]  such 
that 

C 

C'2(to,T')||u(to)||  <  ^o,ki  G  [k,oo). 

Now  we  can  also  select  ko  and  ki  so  large  that 

Thus  we  have 

||u(t)||^  <  <5, 

and  it  follows  that 

lim  ||w(t)l|  =  0  for  any  t  >  0 

ko,ki-^oo 
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and  the  convergence  is  uniform  on  any  interval  [to,T]  with  to  >  0. 

Furthermore,  due  to  statement  6.  above,  we  see  that  the  convergence  takes  place  on  the 
whole  interval  [0,  T]. 

As  an  example  consider  the  controlled  viscous  Burgers’  system  (1.29)  with  disturbance 
/  =  0  and  initial  condition  w{x,  0)  =  — 


t-axis  ®  ®  X-axis  ^  ^  x-axIs 


ko  =  h  =  100  ko  =  ki  =  1000 

The  first  figure  depicts  the  open  loop  system  dynamics,  i.e.,  fco  =  0,  fci  =  0  and  the 
trajectories  of  the  zero  dynamics  system  corresponding  to  fco  =  fci  =  oo.  Note  that  the 
solution  of  the  uncontrolled  open  loop  system  approaches  the  nonzero  steady  state 

w{x,t)  =  m33^{<f>,l)^  [  x'^il-xYdx. 

Jo 

In  general  for  the  Burgers’  system  (1-29)  the  nonzero  equilibrium  guaranteed  by  the  Center 
Manifold  Theorem  appears  to  depend  on  the  value  of  e  unlike  the  linearization  of  the  Burgers’ 
system  (the  heat  equation  with  homogeneous  Neumann  boundary  conditions)  for  which  it 
known  that  the  steady  state  is  always  equal  to  the  integral  of  the  initial  data,  for  all  e.  In 
the  second  and  third  figures,  we  have  depicted  the  dynamics  of  the  closed  loop  system  for 
various  values  of  the  gain  parameters  k  =  ko  =  ki  =  1,10, 100, 1000.  We  have  computed  the 
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mayirmim  error  between  the  trajectories  of  the  closed  loop  dynamics  and  the  zero  dynamics 
on  a  uniform  grid  in  x  —  t  space. 

k  1  10  100  1000 

error  .025  .0088  .0012  .0001 

In  the  next  example,  we  consider  the  controlled  viscous  Burgers’  system  (1.29)  with 
f{x)  =  (0.5)  sin(27ra;)  and  initial  condition  2a;^(l  —  x)^. 
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ko  =  ki  =  100  ko  =  ki  =  1000 


k  1  10  100  1000 

error  .0687  .0220  .0029  .0003 

As  a  final  example  we  consider  a  Burgers’  system  with  distributional  forcing  term  /  (x)  = 
(0.25)  ^3:0(3:)  with  Xo  =  0.75  and  zero  initial  condition. 


ko  =  ki  =  1  ko  =  ki  =  10 
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r 


ko  =  h  =  100  ko  =  ki  =  1000 


k  1  10  100  1000 

error  .144  .018  .002  .0002 

In  the  paper  [91]  we  have  extended  these  results  on  convergence  of  trajectories,  at  least  for 
small  intial  data,  to  a  general  class  of  one  dimensional  convection  reaction  diffusion  equations 
and  we  are  currently  complete  in  the  process  of  extending  these  results,  for  all  finite  energy 
initial  data,  to  a  general  class  of  convective  reaction  diffusion  equations  on  bounded  domains 
in  R”. 


1.3  Boundary  Stabilization  for  Nonlinear  Distributed  Parameter 
Systems 

Motivated  in  part  by  problems  of  flow  control  and  combustion  control,  where  nonlinear  effects 
actually  can  improve  mixing,  there  has  been  considerable  attention  given  in  the  literature 
to  the  study  of  asymptotic  or  steady  state  properties  of  solutions  of  nonlinear  distributed 
parameter  systems,  such  as  Navier-Stokes  equations,  which  contain  both  nonlinear  convec¬ 
tive  terms  and  diffusive  terms.  In  addition,  the  need  to  control  systems  in  which  the  rigid 
aerodynamics  are  coupled  with  the  effects  of  fluid  flow  and  flexible  structural  modes  also  un¬ 
derscores  the  importance  of  developing  a  feedback  design  methodology  capable  of  producing 
controllers  for  nonlinear  distributed  parameter  systems  in  a  systematic  way. 

Of  considerable  interest  then  is  the  question  of  the  effect  of  boundary  control,  modeling 
the  use  of  point  actuators  and  sensors,  in  controlling  or  influencing  the  steady-state  response 
of  forced  nonlinear  distributed  parameter  systems.  In  the  recent  paper  [83]  we  consider  an 
example  of  such  a  boundary  control  problem  for  the  forced  Burgers’  equation  on  a  finite 
interval.  This  remarkable  system  (see  e.g.  [76])  has  been  studied  for  some  time  [77]  and 
was  extensively  developed  by  Burgers  [78]  as  a  simplified  fluid  flow  model  which  nonetheless 


27 


exhibits  some  of  the  important  aspects  of  turbulence.  It  was  later  derived  by  Lighthill  [79] 
as  a  second  order  approximation  to  the  one  dimensional  unsteady  Navier-Stokes  equation. 

The  uncontrolled  problem  studied  in  [83]  corresponds  to  Burgers’  equation  with  homo¬ 
geneous  Neumann  boundary  conditions.  Since  constants  are  global  solutions  we  see  that  the 
uncontrolled  system  is  not  asymptotically  stable.  For  this  reason  the  question  of  feedback 
regularization  becomes  an  important  issue.  In  particular,  we  are  interested  in  introducing 
a  boundary  feedback  control,  delivering  a  closed  loop  system,  for  which  there  is  a  globally 
defined  dynamics  in  L^.  Motivated,  in  part,  by  classical  feedback  design  and  root  locus  for 
finite  dimensional  control  systems,  we  formally  employ  a  simple  proportional  error  feedback 
control  law.  For  the  linerarization  about  zero  of  the  Burgers’  system,  i.e.,  the  one  dimen¬ 
sional  heat  equation,  such  a  control  law  can  be  rigorously  developed  based  on  the  work  in 
[94] .  Indeed,  for  the  linear  problem  this  proportional  error  feedback  exponentially  stabilizes 
the  closed  loop  system  in  all  of  L^.  The  boundary  feedback  control  is  equivalent  to  introduc¬ 
ing  radiation  boundary  conditions  at  each  end  of  the  rod  with  parameters  in  the  boundary 
conditions  considered  as  “gains”  or  control  parameters.  In  our  work  we  have  demonstrated 
that  this  same  feedback  scheme  provides  a  closed  Burgers’  system  with  a  globally  defined 
dynamics  which  for  the  unforced  problem  is  globally  Lyapunov  stable.  Thus  this  example 
demonstrates  that  boundary  control  can  significantly  regularize  nonlinear  distributed  param¬ 
eter  systems;  including  systems  of  hydrodynamic  type  which  contain  nonlinear  convective 
terms. 

For  small  initial  data  in  restrictions  on  the  forcing  terms  can  be  considerably  relaxed 
and  we  show  that  there  is  a  globally  (in  time)  defined  dynamics  possessing  a  local  absorbing 
ball.  For  stationary  forcing  terms  we  have  shown  that  the  corresponding  nonlinear  semi¬ 
groups  define  compact  operators  for  all  positive  time,  from  which  it  was  possible  to  deduce 
the  existence  of  a  local  attractor  with  finite  fractal  and  finite  Hausdorf  dimension.  As  a 
corollary  to  this  analysis,  we  also  showed  that  the  unforced  system  is  locally  exponentially 
stable  in  ,  representing  a  marked  improvement  over  a  similar  result  for  small  initial  data 
in  H^{0, 1)  reported  in  [97]. 

Concerning  the  method  of  proof,  it  is  well-known  [80]  that  the  one-dimensional  Burgers’ 
equation  with  an  external  forcing  term  can  be  reduced  by  the  Hopf-Cole  substitution  to 
the  one-dimensional  heat  equation  with  a  “potential”  term,  which  might  seem  to  offer  a 
straightforward  method  for  asymptotic  analysis.  However,  for  Neumann  or  radiation  bound¬ 
ary  conditions,  this  substitution  introduces  quadratic  nonlinearities  at  the  boundary  of  the 
interval,  in  sharp  contrast  with  the  case  of  Dirichlet  boundary  conditions.  Also,  a  complicat¬ 
ing  feature,  which  arises  for  boundary  conditions  which  are  neither  Dirichlet  nor  periodic, 
is  a  nonlinear  term  which  persists  in  the  energy  balance  relation  and  presents  considerable 
difficulty  in  obtaining  a  priori  estimates  from  which  the  existence  of  an  absorbing  ball  and 
the  convergence  of  Galerkin  approximations  can  be  established  as  in  [73,  81]. 

Our  method  can  be  briefly  characterized  in  the  following  way.  The  state  space  for  our 
system  is  L^(0, 1)  and  for  arbitrary  initial  data  in  1/^(0, 1)  we  first  prove  local  in  time  existence 
and  uniqueness  of  the  weak  solution.  On  this  step  we  use  the  energy  method  and  Galerkin 
approximations  following  closely  the  classical  approach  developed  in  [73,  75].  The  energy 
estimates  obtained  for  general  initial  data  in  this  step  are  local  in  time,  and  for  this  reason 
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we  cannot  prove  global  existence  of  solutions  in  the  first  step.  We  also  stress  that  since  our 
state  space  is  L^(0, 1)  we  cannot  use  well  known  existence  results  for  quasilinear  parabolic 
equations  in  the  spaces  of  Holder  continuous  functions  based  on  the  Leray-Schauder  fixed 
point  theorem  [75]. 

Our  second  step,  which  is  technically  the  most  involved,  is  the  systematic  investigation 
of  the  regularity  properties  of  the  local  solution.  In  other  words,  we  study  the  smoothing 
properties  of  the  local  in  time  dynamics.  We  show,  in  particular,  that  under  appropriate 
assumptions  on  the  forcing  term,  the  weak  solution  is,  in  fact,  classical  for  positive  values  of 
time.  Our  method  also  provides  upper  bounds  for  the  rate  at  which  the  derivatives  of  the 
solution  can  blow  up  as  t  >  0.  Certainly,  the  smoothing  properties  of  parabolic  equations  is 
a  well  known  phenomenon  which  has  been  studied  extensively  in  the  literature.  In  particular, 
the  smoothing  property  for  Burgers’  equation  with  periodic  boundary  conditions  is  very  well 
explained  in  [67].  However,  as  we  show  below,  the  proof  of  this  property  and  derivation  of 
the  existence  of  higher  order  derivatives  of  the  solution  requires  a  significant  effort  when  a 
boundary  feedback  control  is  implemented. 

Finally,  we  combine  the  local  in  time  existence  and  regularity  result  with  a  maximum 
principle  argument  to  obtain  the  global  in  time  existence  and  regularity  of  solutions. 

Yet  another  result  obtained  during  this  research  effort  concerns  the  effect  of  boundary 
control  on  the  structure  of  attractors  for  Burgers’  equation.  In  the,  soon  to  be  published, 
work  [92],  we  prove  that  the  structure  of  the  attractor  is,  in  general,  nontrivial.  Indeed, 
for  a  special  class  of  forcing  terms  we  are  able  to  show  that  as  the  gain  parameters  in  the 
boundary  feedback  control  are  increased  from  small  positive  values  to  large  positive  values, 
the  number  of  stationary  solutions  vary  fi:om  three  to  one.  This  is  in  keeping  with  our  results 
on  conergence  of  attractors,  discussed  in  the  last  subsection,  where  is  can  be  shown  that  the 
zero  dynamics  systems  has  a  single,  global,  asymptotically  stable,  equilibrium. 

1.4  Output  Regulation  for  Linear  Distributed  Parameter  Systems 

In  general,  the  ability  to  systematically  control  or  influence  nonlinear  effects  would  make 
a  substantial  contribution  to  existing  and  emerging  commercial  and  defense  research  and 
development  programs.  Notable  examples,  widely  appreciated  within  the  aerospace  industry, 
include  the  development  of  flight  controllers  for  high  angle-of-attack  or  high  agility  aircraft. 
Indeed,  the  importance  of  including  the  nonlinear  behavior  of  aerodynamic  parameters,  such 
as  the  coefficient  of  lift,  as  a  function  of  the  angle-of-attack  has  long  been  recognized  since 
at  high  angles-of-attack,  wind  angle  moments  also  exhibit  nonlinear  effects  which  cannot  be 
ignored.  Another  area  of  interest  is  the  control  of  flutter,  which  can  shorten  the  life  cycle  of 
aircraft  and  aircraft  parts.  Indeed,  one  example  of  the  potential  impact  of  nonlinear  control 
in  problems  of  flow  control  is  in  the  control  of  instabilities  in  the  unsteady  separated  shear 
layer,  which  has  been  experimentally  shown  to  greatly  influence  stall  and  lift  behavior  at 
high  angles  of  attack. 

It  is  worth  noting,  however,  that  the  active  control  in  experiments,  such  as  Batill  and 
Mueller  [69],  is  based  on  a  priori  harmonic  forcing,  while  in  nonlinear  systems  with  resonance 
it  is  known  that  simple  harmonic  forcing  will  not  necessarily  produce  the  desired  response 
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(see  e.g.  [12]  for  an  analysis  of  the  steady  state  response  of  nonlinear  systems  by  center 
manifold  methods).  This  point  was  also  illustrated  theoretically  by  Keefe  [30],  who  showed 
that  the  success  of  a  priori  control  for  the  Ginzburg- Landau  equation  was  dependent  on  the 
initial  state  of  the  system,  which  of  course  may  not  be  controlled  or  even  known,  and  that 
therefore  undesirable  responses  are  to  be  expected  in  the  nonlinear  regime.  Moreover,  the 
computation  by  Fuglsang  and  Cain  [9]  of  flow  over  an  open  cavity  suggests  that  harmonic 
forcing  at  non-resonant  harmonic  frequencies  can  produce  a  limit  cycle  or  chaotic  response 
that  is  far  more  severe  than  the  natural  harmonic  resonance. 

One  of  the  central  problems  in  control  theory  is  to  control  a  flxed  plant  in  order  to  have 
its  output  track  a  reference  signal  (and/or  reject  a  disturbance)  produced  by  an  external 
generator  or  exogenous  system.  Generally  two  versions  of  this  problem  are  considered.  In  the 
first,  the  state  feedback  regulator  problem,  the  controller  is  provided  with  full  information 
of  the  state  of  the  plant  and  exosystem.  For  the  second,  more  realistic  error  feedback 
regulator  problem,  only  the  components  of  the  error  are  available  for  measurement.  For  linear 
finite  dimensional  systems  it  has  been  shown  by  Francis  [104]  that  the  solvability  of  the  the 
regulator  problems  is  equivalent  to  the  solvability  of  a  pair  of  linear  matrix  equations.  This 
in  turn  can  be  characterized  as  a  property  of  the  transmission  polynomials  of  the  composite 
system  formed  from  the  plant  and  the  exosystem  as  it  was  shown  by  Hautus  [106].  Francis 
and  Wonham  [105]  have  also  shown  that  any  regulator  that  solves  the  error  feedback  problem 
has  to  incorporate  a  model  of  the  exogenous  system  generating  the  reference  signal  which  is 
to  be  tracked  and/or  the  disturbance  that  must  be  rejected.  This  property  is  known  as  the 
internal  model  principle. 

Similar  results  have  been  established  for  finite  dimensional  nonlinear  systems  in  [12]  (see 
also  [107])  in  case  the  plant  is  exponentially  stabilizable  and  the  exosystem  has  bounded 
trajectories  that  do  not  trivially  converge  to  zero.  In  particular,  it  is  shown  in  [12]  that 
the  solvability  conditions  given  by  Francis  in  the  linear  case  can  be  naturally  generalized  to 
the  solvability  of  a  pair  of  nonlinear  equations  -  the  regulator  equations  —  that  express  the 
existence  of  a  local  manifold  on  which  the  actual  and  reference  outputs  coincide  and  which 
can  be  rendered  invariant  using  feedback.  The  paper  [12]  also  contains  an  interpretation 
Hautus’  result  in  the  nonlinear  setting  in  terms  of  the  notion  of  zero  dynamics.  They 
expressed  the  solvability  conditions  as  a  special  property  of  the  zero  dynamics  of  the  system 
composed  from  the  plant  and  the  exosystem  with  the  tracking  error  as  its  output.  In  [12]  it 
is  proven  that  the  nonlinear  regulator  equations  are  solvable  if  and  only  if  the  zero  dynamics 
of  this  composite  system  contains  a  controlled  invariant  submanifold  on  which  the  dynamics 
is  a  diffeomorphic  copy  of  the  flow  of  the  exosystem. 

In  our  forthcoming  work  [102]  we  have  developed  a  geometric  theory  of  output  feedback 
regulation  for  infinite  dimensional  linear  control  systems  assuming  that  the  control  and 
observation  operators  are  bounded.  In  particular  we  have  obtained  results  characterizing 
the  solvability  of  both  state  and  error  feedback  regulator  problems  in  terms  of  the  solvability 
of  certain  equations  referred  to  as  the  regulator  equations.  The  main  difficulties  that  arise 
in  extending  the  work  in  [12]  to  the  distributed  parameter  case  are  obvious:  the  phase 
space  is  infinite  dimensional;  the  state  operator  is  unbounded  and  consequently  only  densely 
defined;  there  is  no  direct  analog  of  the  Jordan  decomposition  consequently  some  care  must 
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be  exercised  in  dealing  with  the  spectra  of  certain  composite  systems. 

To  motivate  our  development  of  the  linear  regulator  theory  for  bounded  inputs  and 
outputs,  let  us  consider  a  simple  example  in  which  we  have  a  plant  consisting  of  a  one 
dimensional  heat  equation  on  a  finite  rod  with  end  points  insulated,  i.e.,  Neumann  boundary 
conditions.  Our  specific  objective  is  to  design  a  feedback  law  which,  for  arbitrary  initial  data, 
will  drive  the  output  -  consisting  of  the  average  temperature  about  a  fixed  point  on  the  rod 
-  in  order  to  track  a  periodic  temperature  profile. 

For  this  example,  the  first  of  our  three  basic  ingredients,  the  existence  of  a  steady  state 
response  is  classical  and  for  the  second  ingredient,  stabilizability,  the  plant  is  easily  stabilized 
by  a  low  order  feedback.  For  this  example,  the  third  ingredient,  shaping  the  response,  is  most 
difiicult.  Once  we  have  stabilized  the  system,  driving  the  system  with  a  sinusoid  produces  an 
output  which  converges  rapidly  to  a  periodic  motion,  but  unfortunately,  the  amplitude  and 
phase  are  not  correct.  The  most  important  design  feature  is  shaping  the  output  by  choosing 
an  appropriate  periodic  input.  However,  applying  the  regulator  theory  discussed  below,  the 
tracking  problem  is  solved  immediately. 

Consider  a  controlled  SISO  one  dimensional  heat  equation  on  the  interval  [0, 1]  : 


(1.1) 

||(0,.)-0.  |(U)  =  0 

(1.2) 

z{x,  0)  =  (j){x), 

y{t)  =  Cz. 

(1.3) 

We  can  formulate  the  system  (1.1),  (1.3)  in  the  usual  {A,B,C)  form  by  choosing 

z{t)  =  Az{t)  +  Bu{t), 

y{t)  =  Cz{t),  (measured  output) 


2(0)  =  Zo, 

where 

1.  The  operator  A  =  cP/dx^  with  Neumann  boundary  conditions  defines  an  unbounded 
selfadjoint  operator  in  Z  specified  by  the  dense  domain, 

V{A)  =  {ct>e  H\0, 1)  :  <^'(0)  =  =  0}. 

2.  The  input  corresponds  to  a  spatially  uniform  temperature  input  over  a  small  interval 
about  a  fixed  point  xq  €  (0, 1): 

Bu  =  h{x)u  (1.4) 
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3.  The  output  corresponds  to  the  average  temperature  over  a  small  interval  about  a  point 

Xi  G  (0, 1). 

C(l)=  (  c(x)(j){x)dx  (1.5) 

Jo 

c{x)  =  ■  (1-6) 

and  where  Z  =  L^(0, 1)  is  the  state  of  the  system,  a  separable  Hilbert  space  (state  space) 
with  inner  product 

<^,ip>=  /  ^{x)'4){x)dx. 

Jo 

The  spectrum  of  A  consists  of  pure  point  spectrum 

a{A)  =  {-k'^TT^jZo 

with  corresponding  orthonormal  eigenvectors 

i^oix)  =  1  and  '4’kix)  =  V2cos(knx)  for  A:  >  1. 

In  this  example,  we  consider  the  problem  of  tracking  a  sinusoid,  yr{t)  =  M  sin(Q:t).  Thus 
for  a  suitable  initial  condition  yr{t)  is  the  output  of  a  harmonic  oscillator: 


W  : 

=  Sw, 

0 

a 

“1 

—a 

:  0 

t(;(0) 

= 

1  i 

yr  = 

[1 

Ojw 

The  ultimate  design  objective  for  this  problem  would  be  to  find  a  control  u  that  will  force 
the  output  y{t)  to  track  the  reference  trajectory  yr,  i.e.,  so  that  the  error  satisfies 

e{t)  =  y(t)  —  yr{t)  0  ajs  t  ^  oo. 


In  order  to  demonstrate  the  importance  of  stabilization  and  shaping  the  response  we  first 
consider  choosing  an  input  u  =  yr  =  Msm{at),  i.e.,  driving  the  system  with  the  desired 
harmonic  reference  signal.  As  might  be  expected,  due  to  the  zero  eigenvalue  which  causes 
the  original  uncontrolled  system  not  to  be  asymptotically  stable,  the  controlled  output  does 
not  oscillate  about  zero. 

For  our  numerical  simulation,  we  have  chosen  Xq  =  3/4,  =  1/4,  uq  =  ui  =  1/4, 

M  =  1,  a  =  2  with  initial  condition  4){x)  =  4a;^(3/2  —  a;)  and  we  have  set  the  control  input 
u  —  sin(2t).  Figure  1  contains  a  plot  of  the  resulting  outputs  y  and  yr  while  Figure  2  depicts 
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the  error  e  —  y  —  Ur  and  Figure  3  depicts  the  entire  solution  surface  for  the  heat  problem 
due  to  this  simple  harmonic  forcing. 


Fig.  1:  y  (solid  line),  yr  (dashed  line) 


Fig.  2:  eit)  for  harmonic  forcing 
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Next  we  incorporate  a  stabilizing  feedback  term  in  the  control.  Thus  we  consider  a 
feedback  control  in  the  form  u  =  Kz  +  Msin(2t)  where  u  =  Kz  is  a,  stabilizing  feedback.  A 
rank  one,  bounded  stabilizing  feedback  operator  for  this  problem  is  given  by  (cf.  [103]) 

Kcj)  =  —(3  <  (f>,l  >,  /3  >  0. 

Indeed  it  can  be  shown  that  {A  +  BK)  is  a  discrete  Reisz  spectral  operator  with  simple 
spectrum  given  by  {— /?}  U  Iii  Figures  4  -  6  we  have  again  plotted  the  resulting 

outputs,  error  and  solution  surface,  respectively,  with  =  0.2. 


Fig.  4:  y  (solid  line),  yr  (dashed  line) 
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Fig.  5:  e(t)  harmonic  forcing  with  stability 


The  stabilizing  feedback  has  provided  an  output  which  now  appears  to  converge  to  a 
periodic  trajectory  about  zero,  as  desired,  but  the  resulting  asymptotic  amplitude  and  phase 
are  not  correct. 

Therefore,  as  suggested  in  the  statement  of  the  state  feedback  regulator  problem  (given 
in  the  next  section)  we  consider  finding  a  feedback  law  in  the  more  general  form 

u  =  Kz  +  Lw 

where  K  is  a,  stabilizing  feedback  operator  and  L  =  (F  —  KT[)  with  11  =  [Hi,  112]  £  £(]R^,  Z) 
and  r  =  [71,72]  £  IR^  are  solutions  of  the  so-called  “regulator  equations”  (described  in 
Theorem  1  below.)  In  this  example  the  first  regulator  equation  reduces  to  a  pair  of  coupled 
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second  order  ordinary  differential  equations  with  boundary  conditions, 


i’ 


III  (a;)  +  ocIl2{x)  =  -b{x)'yi  (1.7) 

112(0:)  -  Q;ni(a;)  =  —h{x)^2- 

n;(o)  =  n;(i)  =  o 

n'2(o)  =  n' (1)  =  0. 

The  boundary  conditions  in  (1.7)  correspond  to  the  requirements  that  Hi  and  112  he  in  the 
domain  of  A. 

The  parameters  71  and  72  are  chosen  to  satisfy  a  second  regulator  equation,  which  in  this 
case  reduces  to  the  additional  constraints 


<c,ni>=l,  <c,n2>=0.  (1.8) 

These  equations  can  be  solved  off-line  either  analytically  or  numerically.  The  resulting 
outputs  and  solution  surface  for  the  system  with  control  u  =  Kz+Lw  are  depicted  in  Figures 
7  and  8.  In  Figure  9  we  have  plotted  the  error  e{t)  =  y{t)  —  yr{p)  and  note  that  the  error 
appears  to  tends  to  zero.  Indeed,  as  a  result  of  the  center  manifold  methods  applied  in  our 
proofs,  we  expect  that  the  error  tends  to  zero  exponentially. 
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0  0 

Fig.  8:  Solution  Surface  for  u  =  Kz  +  Lw 


Fig.  9:  e{t)  using  regulator  equations 


Regulator  Problem 

Consider  a  plant  described  by  an  abstract  distributed  parameter  control  system  in  Hilbert 
space: 


z(t)  =  Az{t)  +  Bu{t)  +  d{t),  (1.9) 

=  Cz{t),  (measured  output)  (1-10) 

2(0)  =  Zo, 

where  z  E  Z  is  the  state  of  the  system,  Z  is  a  separable  Hilbert  space  (state  space),  u  e  U 
is  an  input  ,  y  E  Y  is  the  measmed  output,  U  and  Y  are  Hilbert  control  and  output 
spaces,  respectively.  A  is  assumed  to  be  the  infinitesimal  generator  of  a  strongly  continuous 
semigroup  on  a  Hilbert  space  Z,  B  E  C{JJ,  Z)  and  C  E  C{Z,Y).  Here  we  use  the  notation 
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C{Wi^W2)  to  denote  the  set  of  all  bounded  linear  operators  from  a  Hilbert  space  Wi  to  a 
Hilbert  space  W2-  The  term  d{t)  represents  a  disturbance. 

In  addition  we  will  assume  that  there  exists  a  finite  dimensional  linear  system,  referred 
to  as  the  exogenous  system  (or  exosystem),  that  produces  a  reference  output  yr{t)  and  which 
is  also  used  to  model  the  disturbance  d[t)\ 


w{t)  =  Sw{t) 

(1.11) 

yr{t)  =  Qw{t) 

(1.12) 

d{t)  =  Pw{t) 

(1.13) 

O 

II 

(1.14) 

Here  S  G  £(E^),  Q  G  £(R^  Y)  and  P  G  £(R^  Z). 

We  will  refer  to  the  difference  between  the  measured  and  reference  outputs  as  the  error 

e(t)  =  y{t)  -  yr{t)  =  Cz{t)  -  Qw{t).  (1.15) 


I.  Linear  State  Feedback  Regulator  Problem: 


exo-system 


z 


STATE  FEEDBACK 

Find  a  feedback  control  law  in  the  form  u{t)  —  Kz{t)  +  Lw{t)  such  that  K  G  JC{Z,  U), 
L  G  £(R^  U)  and 

I.a  z{t)  =  {A  +  BK)z{t)  is  stable,  i.e.  {A  +  BK)  is  the  infinitesimal  generator  of  an 
exponentially  stable  Co  semigroup. 

I.b  For  the  closed  loop  system 

z{t)  =  (^  +  BK)z{t)  +  {BL  +  P)w{t), 

w{t)  =  Sw{t), 
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the  error 
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e(t)  —  Cz{t)  —  Qw{t)  — »•  0  as  t  ^  oo, 
for  any  initial  condition  in  Z  x 


II.  Linear  Error  Feedback  Regulator  Problem: 

exo-system 


ERROR  FEEDBACK 

Find  an  error  feedback  controller  of  the  form 

X{t)  =  FX{t)  +  Ge{t), 
u{t)  =  HX{t) 

where  X{t)  G  X  for  t  >  0,  X  is  a  Hilbert  space,  G  G  C{Y,X),  H  G  C{X,U)  and  F  is  the 
infinitesimal  generator  of  a  Go  semigroup  on  X  with  the  properties  that 

Il.a  The  system 


z{t)  =  Az{t)  +  BHX{t), 
X{t)  =  FX{t)  +  Ge(t) 


is  exponentially  stable  when  w  =  0,  i.e. 
an  exponentially  stable  Gq  semigroup,  and 

II. b  for  the  closed  loop  system 


A  BH 
GC  F 


is  the  infinitesimal  generator  of 


z{t)  =  Az{t)  +  BHX{t)  +  Pw{t), 
X(t)  =  GGz(t)  +  FX(t)  -  GQw(t), 


w(t)  —  Sw{t) 
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the  error 

e{t)  =  Cz{t)  —  Qw{t)  — >  0  ast^oo, 
for  any  initial  condition  in  Z  x  X  xEf. 

We  impose  the  following  standard  assumptions. 

Three  basic  assumptions: 

HI  Without  loss  of  generality,  we  assume  that  cr(5)  is  contained  in  the  closed  right  half 
plane.  (Modes  corresponding  to  spectrum  in  the  open  left  half  plane  decay  exponen¬ 
tially  to  zero.) 

H2  The  pair  {A,  B)  is  exponentially  stabilizable,  i.e  there  exists  K  €  C{Z,  U)  such  that 
A  BK  is  the  infinitesimal  generator  of  an  exponentially  stable  Co  semigroup, 

H3  The  pair 


is  exponentially  detectable,  i.e  there  exists  G  €  C{Y,  Z  x  R^), 

C=  [  1 ,  G,eC{Y,Z),  C2e£(y,R'') 

such  that 


is  the  infinitesimal  generator  of  an  exponentially  stable  Co  semigroup. 

The  first  main  results  from  [102]  characterizing  the  solvability  of  the  regulator  problems 
for  linear  distributed  parameter  systems  are  given  in  Theorems  1  and  2,  respectively. 

THEOREM  1:  Let  HI  and  H2  hold.  The  linear  state  feedback  regulator  problem  is  solvable 
if  and  only  if  there  exist  mappings  11  G  Z)  with  RanfR)  C  T*(A)  and  T  €  U) 

satisfying  the  “regulator  equations,  ” 

ns  =  Au  +  Br  +  p,  (1.16) 

Cn  =  Q.  (1.17) 

In  this  case  a  feedback  law  solving  the  state  feedback  regulator  problem  is  given  by  u{t)  = 
Kz{t)  -h  (r  —  KlT)w{t),  where  K  is  any  exponentially  stabilizing  feedback  for  {A,  B). 
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THEOREM  2:  Let  HI,  H2  and  H3  hold.  The  linear  error  feedback  regulator  problem  is 
solvable  if  and  only  if  there  exist  mappings  11  G  Z)  and  F  G  U)  with  i?an(n)  C 

V{A),  such  that 

US  =  AU  +  Br  +  P  (1.18) 


cn  =  Q.  (1.19) 

With  this  n  and  F  a  controller  solving  the  error  feedback  regulator  problem  is  given  by 

X{t)  =  FX{t)  +  Ge{t), 
u(t)  =  HX{t). 


where  X  £  X  =  Z  y. 


where 


G  = 


Gi 

G2 


,  H=  [K  (F-l^n)], 


F  = 


Fn  Fi2 
F21  F22 


Fn  =  {A  +  BK-GiC), 

Fi2  =  {P  +  5(r  -  KU)  +  GiQ) , 

F21  =  —G2C, 

F22  =  (<5  +  G2Q) 

and  K  G  C{Z,  U)  is  an  exponentially  stabilizing  feedback  for  the  pair  (A,  B)  and 
an  exponentially  stabilizing  output  injection  (such  K  and  G  exist  by  H2  and  H3j. 

Solvability  Criteria 

We  have  the  following  solvability  conditions  for  the  regulator  equations  for  SISO  systems 
which  relates  solvability  to  the  transmission  zeros  of  the  system  (cf.  [101]). 


G2 
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THEOREM  3:  For  a  SISO  system,  if  the  exo-system  is  neutrally  stable,  the  regulator 
equations  are  solvable  if  and  only  if  the  nonresonance  condition 

g{\)  0  for  all  A  G  a{S) 


is  satisfied.  Here 


g{s)  =  C{sI-A)-^B 


is  the  transfer  function  of  the  plant.  That  is,  if  and  only  if,  the  eigenvalues  of  S  are  not 
poles  or  transmission  zeros  of  the  plant. 
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As  a  simple  example  of  this  result,  let  us  return  to  our  controlled  heat  equation  for  which 
the  transfer  function  can  be  written  as 

g{s)  =  c{sl  —  A)~^b,  s  E  C. 


In  our  specific  numerical  example 

_  2sinh(-v/i/2) 

s-yicosh(-ys/2) ' 

More  generally,  defining  the  operators 

Ri  =  + 


i?2  =  Q:(A^  +  Q!^)  \ 


we  can  explicitly  express  the  solutions  IIi,  112,  7i  and  72  of  the  regulator  equations  in  (1.7) 

as  T.,  /.  N 

Reg{ia) 

“  \g{ia)\^  ’ 


72  = 


Im^(iQ:) 


and 


111  I  /  •  1  _  M9  “^2^5 


\g{iaW 


Reg{ia)  lmg{ia) 

^^2  =  I  /■  M9'^20  +  I  /•  MO  -^iP- 

For  a  ^  0  the  conditions  of  Theorem  3  are  satisfied  and  all  the  expressions  make  sense. 
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Author:  W.  Lin  and  C.I.  Byrnes 

Publisher:  Elsevier  Science  Ltd 

Volume:  31,  Page:  269-274,  Year  Published:  1995. 

Name  of  Journal:  Automatica 

Title  of  Article:  “Discrete-Time  Nonlinear  H-infinity  Control  with  Measurement  Feed¬ 
back” 

Author:  W.  Lin  and  C.I.  Byrnes 

Publisher:  Elsevier  Science  Ltd 

Volume:  31,  Page:  419-434,  Year  Published:  1995. 

Name  of  Journal:  Systems  and  Control  Letters 

Title  of  Article:  “Remarks  on  Linearization  of  Discrete-Time  Autonomous  Systems 
and  Nonlinear  Observer  Design” 

Author:  W.  Lin  and  C.I.  B5rrnes 

Publisher:  Elsevier  Science  BV 

Volume:  25,  Page:  31-40,  Year  Published:  1995. 

Name  of  Journal:  IEEE  Trans.  Automat.  Contr. 

Title  of  Article:  “H-infinity  Control  of  Discrete-Time  Nonlinear  Systems” 

Author:  W.  Lin  and  C.I.  Byrnes 
Publisher:  IEEE 

Volume:  41,  Page:  494-510,  Year  Published:  1996. 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 

Title  of  Article:  “  H-infinity  Control  of  Nonlinear  Systems  with  Sampled  Measurements” 

Author(s):  S. Suzuki,  A.Isidori  and  T.J.Tarn 

Publisher:  Birkhauser 

Volume:  5  ,  Page(s):  259-262, 

Month  Published:  March,  Year  Published:  1995 

Name  of  Journal:  IEEE  Transactions  on  Automatic  Control 

Title  of  Article:  “H-infinity  control  of  discrete  time  nonlinear  control  systems  via  state 
and  full  information  feedback” 

Author:  C.I.Byrnes  and  W.  Lin 
Publisher:  IEEE  Volume:  ,  Page: 

Year  Published:  1996  (to  appear). 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 
Title  of  Article:  “Solution  of  Nonlinear  Lagrange  and  Bolza  problems  via  Riccati  par- 
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tial  differential  equations” 

Author:  C.I. Byrnes 

Publisher:  ,  Page:  ,  Year  Published:  . 

Name  of  Journal:  IEEE  Transactions  on  Automatic  Control 

Title  of  Article:  “Robust  regulation  for  nonlinear  systems  with  gain-bounded  uncer¬ 
tainties” 

Author:  A.Isidori  and  T.J.Tarn 
Publisher:  IEEE 

Volume:  40,  Page:  1744-1754,  Year  Published:  1995. 

Name  of  Journal:  Automatica 

Title  of  Article:  “Structurally  Stable  Output  Regulation  of  Nonlinear  Systems” 

Author:  C.I. Byrnes,  F.  Belli  Priscoli,  A.  Isidori  and  W.Kang 
Publisher:  Pergamon  Press 

Volume:  no.  3,  Page:369-385,  Year  Published:  1997. 

Name  of  Journal:  IEEE  Control  Systems  Magazine 

Title  of  Article:  “Force  regulation  and  contact  transition  control” 

Author:  T.J.Tarn,  Y.  Wu,  N.Xi  and  A.Isidori 
Publisher:  IEEE 

Volume:  16,  Page:  32-40,  Year  Published:  1996. 

Name  of  Journal:  Systems  and  Control  Letters 

Title  of  Article:  “A  note  on  almost  disturbance  decoupling  for  nonlinear  minimum 
phase  systems” 

Author:  A.Isidori 

Publisher:  Elsevier  Science  BV  Volume:  27,  Page:  191-194, 

Year  Published:  1996. 

Name  of  Journal:  IEEE  Trans.  Automat.  Contr. 

Title  of  Article:  “H-infinity  Control  of  Discrete-Time  Nonlinear  Systems” 

Author:  W.  Lin  and  C.  I.  Byrnes 
Publisher:  IEEE 

Volume:  40,  Page:  494-510,  Year  Published:  1996. 

Name  of  Journal:  Automatica 

Title  of  Article:  “Global  Asymptotic  Stabilization  of  General  Nonlinear  Systems  with 
Stable  Rree  Dynamics  via  Passivity  and  Bounded  Feedback” 

Author:  W.  Lin 

Publisher:  Pergamon  Press 

Volume:  32,  Page:  915-924,  Year  Published:  1996. 
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Name  of  Journal:  International  Journal  of  Control 

Title  of  Article:  “Mixed  H2/H-infinity  Control  via  State  Feedback  for  Nonlinear  Sys¬ 
tems” 

Author:  W.  Lin 

Publisher:  Taylor  &  Francis  Ltd 
Volume:  Page:  Year  Published:  1996. 

Name  of  Journal:  Systems  and  Control  Letters 

Title  of  Article:  “Further  results  on  Global  Stabilization  of  Discrete  Nonlinear  Sys¬ 
tems” 

Author:  W.  Lin 

Publisher:  Elsevier  Science  BV 

Volume:  28,  Page:  Year  Published:  1996. 

Name  of  Journal:  Differential  and  Integral  Equations 

Title  of  Article:  “Global  Solvability  for  Damped  Abstract  Nonlinear  Hyperbolic  Sys¬ 
tems” 

Author:  H.T.  Banks,  D.S.  Gilliam,  V.I.  Shubov 
Publisher: 

Volume:  10  no.  2,  Page:  309-332,  Year  Published:  1997. 

Name  ofJournal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 

Title  of  Article:  “High  Gain  Limits  of  Trajectories  and  Attractors  for  a  Boundary 

Controlled  Viscous  Burgers’  Equation” 

Author:  C.I.  Byrnes,  D.S.  Gilliam,  V.I.  Shubov 
Publisher: 

Volume:  ,  Page:  ,  Year  Published:  1996 

Name  of  Journal:  Journal  of  Dynamical  and  Control  Systems 

Title  of  Article:  “On  the  Global  Dynamics  of  a  Controlled  Viscous  Burgers’  Equation” 

Author:  C.I.  Byrnes,  D.S.  Gilliam,  V.I.  Shubov 

Publisher: 

Volume:  ,  Page:  ,  Year  Published:  1996. 

Name  of  Journal:  Automatica 

Title  of  Article:  “Structurally  stable  output  regulation  of  nonlinear  systems” 

Author:  C.I.Byrnes,  F.Delli  Priscoli,  A.  Isidori,  W.Kang 

Publisher:  Elsevier  Science  BV 

Volume:  33 

Page:  369-385 

Year  Published:  1997. 

Name  of  Journal:  IEEE  Trans,  on  Automatic  Control 
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Title  of  Article:  “A  remark  on  the  problem  of  semiglobal  nonlinear  output  regulation” 

Author:  A.  Isidori 

Publisher:  IEEE  Press 

Volume:  43 

Page:  to  appear 

Year  Published:  1997. 

Name  of  Journal:  IEEE  Trans,  on  Automatic  Control 

Title  of  Article:  “On  the  partial  stochastic  realization  problem” 

Author:  C.I.  Byrnes  and  A.  Lindquist 
Publisher:  IEEE  Press 
Volume:  42,  No.  8 
Page:  1049  -  1070 
Year  Published:  1997. 

Name  of  Journal:  SIAM  J.  Cont.  and  Opt. 

Title  of  Article:  “A  convex  optimization  approach  to  the  rational  covariance  extension 
problem” 

Author:  C.I.  Byrnes,  S.V.  Gusev,  and  A.  Lindquist 
Publisher:  SIAM  Press 
Status:  Submitted 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 
Title  of  Article:  “Stationary  solutions  for  a  boundary  controlled  viscous  Burgers’  equa¬ 
tion” 

Author:  A.  Balogh,  D.S.  Gilliam  and  V.I.  Shubov 
Status:  Submitted 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 
Title  of  Article:  “Numerical  stationary  solutions  for  a  viscous  Burgers’  equation” 
Author:  A.  Balogh,  J.A.  Burns,  D.S.  Gilliam  and  V.I.  Shubov 
Publisher:  Birkhauser 
Accepted:  December  1996. 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 
Title  of  Article:  “Harmonic  forcing  for  linear  distributed  parameter  systems” 

Author:  C.I.  Byrnes,  D.S.  Gilliam,  1.  Lauko  and  V.I.  Shubov 
Publisher:  Birkhauser 
Accepted:  December  1996. 

Name  of  Journal:  Journal  of  Mathematical  Systems,  Estimation  and  Control 
Title  of  Article:  “On  the  Riccati  Partial  Differential  Equation  for  Nonlinear  Bolza  and 
Lagrange  Problems” 
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Author:  C.I.  Byrnes 
Publisher:  Birkhauser 
Accepted:  May  1997. 


b.  Peer  Reviewed  Conference  Proceedings: 

Name  of  Conference:  IEEE  Conference  on  Decision  and  Control 
Title  of  Article:  “High  Gain  Limit  for  Boundary  Controlled  Convective  Reaction  Dif¬ 
fusion  Equations” 

Author:  C.I.  Byrnes,  D.S.  Gilliam,  N.  Okasha,  V.I.  Shubov 
Publisher:  IEEE 

Volume:  34th  CDC,  Page:  ,  Year  Published:  1995. 

Name  of  Conference:  IEEE  Conference  on  Decision  and  Control 

Title  of  Article:  “Well-Posedness  for  Controlled  Nonlinear  Damped  Membranes  with 

Fixed  Boimdary” 

Author:  H.T.  Banks,  D.S.  Gilliam,  V.I.  Shubov 
Publisher:  IEEE 

Volume:  34th  CDC,  Page:  ,  Year  Published:  1995. 

Name  of  Conference:  IEEE  Conference  on  Decision  and  Control 

Title  of  Article:  “Global  Robust  Stabilization  of  Minimum-Phase  Nonlinear  Systems” 

Author:  W.  Lin 

Publisher:  IEEE 

Volume:  34th  CDC,  Page:  521-526,  Year  Published:  1995. 

Name  of  Conference:  IEEE  Conference  on  Decision  and  Control 
Title  of  Article:  “Global  Stabilization  of  Discrete  Non- Affine  Systems” 

Author:  W.  Lin 
Publisher:  IEEE 

Volume:  34th  CDC,  Page:  510-511,  Year  Published:  1995. 

Name  of  Conference:  IEEE  Conference  on  Decision  and  Control 
Title  of  Article:  “Mixed  H2/H- infinity  Control  of  Nonlinear  Systems” 

Author:  W.  Lin 
Publisher:  IEEE 

Volume:  34th  CDC,  Page:  333-338,  Year  Published:  1995. 

Name  of  Conference:  13th  IFAC  World  Congress 

Title  of  Article:  “Passivity,  Bounded  State  Feedback  and  Global  Stabilization  of  Non¬ 
linear  Systems” 

Author:  W.  Lin 
Publisher:  Elsevier 

Volume:  E,  Page:  103-108,  Year  Published:  1996. 


55 


Name  of  Conference:  13th  IFAC  World  Congress 

Title  of  Article:  “Solutions  to  the  Output  Regulation  Problem  of  Linear  Singular  Sys¬ 
tems” 

Author:  W.  Lin  and  L.  Dai 
Publisher:  Elsevier 

Volume:  D,  Page:  97-102,  Year  Published:  1996. 

Name  of  Conference:  1995  American  Control  Conference 

Title  of  Article:  “Some  New  Results  on  Stability  and  Observability  of  Discrete-time 
Autonomous  Systems” 

Author(s):  W.  Lin  and  C.I.  Byrnes 
Publisher:  IEEE 

Volume:  6,  Page(s):  4214-4218,  Month  Published:  June, 

Year  Published:  1995 

Name  of  Conference:  IFAC  Nonlinear  Control  Systems  Symposium 
Title  of  Article:  “On  the  Dynamics  of  Boundary  Controlled  Nonlinear  Distributed  Pa¬ 
rameter  Systems” 

Author(s):  C.I.  Byrnes,  D.S.  Gilliam  and  Victor  I.  Shubov  Publisher:  Pergamon 
Volume:  62  Page(s):  913-913,  Month  Published:  June, 

Year  Published:  1995 

Name  of  Conference:  34th  IEEE  Conference  on  Decision  and  Control 

Title  of  Article:  “Well-Posedness  for  Controlled  Nonlinear  Damped  Membranes  with 

Fixed  Boundary  ” 

Author(s):  H.T.  Banks,  D.S.  Gilliam  and  Victor  1.  Shubov 
Publisher:  IEEE 

Volume:  Page(s):  ,  Month  Published:  December, 

Year  Published:  1995 

Name  of  Conference:  35rd  IEEE  Conf.  Decision  and  Control  (Kobe,  Japan,  December 
1996) 

Title  of  Article:  “Global  normal  forms  for  MIMO  nonlinear  systems,  with  application 
to  stabilization  and  disturbance  attenuation” 

Author:  B. Schwartz,  A.  Isidori  and  T.J.  Tarn 
Page:  1041-1046 
Year  Published:  1996. 

Name  of  Conference:  35rd  IEEE  Conf.  Decision  and  Control  (Kobe,  Japan,  December 
1996) 

Title  of  Article:  “L2  disturbance  attenuation  and  performance  bounds  for  linear  non¬ 
minimum  phase  square  invertible  systems” 
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Author:  B. Schwartz,  A.  Isidori  and  T.J.  Taxn 

Page:  227-228 

Year  Published:  1996. 

Name  of  Conference:  4th  European  Control  Conference  (Brussels,  Belgium,  July  1997) 
Title  of  Article:  “Performance  bounds  for  disturbance  attenuation  in  nonlinear  non¬ 
minimum-phase  systems  ” 

Author:  B. Schwartz,  A.  Isidori  and  T.J.  Tarn 
Page:  to  appear 
Year  Published:  1997. 

Name  of  Conference:  36th  IEEE  Conf.  Decision  and  Control  (San  Diego,  CA,  De¬ 
cember  1997) 

Title  of  Article:  “Theoretical  and  numerical  results  for  parameter  estimation  in  non¬ 
linear  elastomers” 

Author:  H.T.  Banks,  D.S.  Gilliam  and  G.  Pinter 
Page:  3733-3738. 

Year  Published:  1997 

Name  of  Conference:  36th  IEEE  Conf.  Decision  and  Control  (San  Diego,  CA,  De¬ 
cember  1997) 

Title  of  Article:  “Output  Regulation  for  Parabolic  Distributed  Parameter  Systems: 
Set  Point  Control” 

Author:  C.I.  Byrnes,  D.S.  Gilliam,  I.  Lauko  and  V.I.  Shubov 
Page:  2231-2236 
Year  Published:  1997 

Name  of  Conference:  36th  IEEE  Conf.  Decision  and  Control  (San  Diego,  CA,  De¬ 
cember  1997) 

Title  of  Article:  “Some  Recent  Results  on  Feedback  Regularization  of  Navier-Stokes 
Equations” 

Author:  A.  Balogh,  D.S.  Gilliam,  and  V.I.  Shubov 
Page:  2231-2236 
Year  Published:  1997 

Name  of  Conference:  36th  IEEE  Conf.  Decision  and  Control  (San  Diego,  CA,  De¬ 
cember  1997) 

Title  of  Article:  “Global  L2-gain  State  Feedback  Design  for  a  Class  of  Nonlinear  Sys¬ 
tems” 

Author:  A.Isidori  and  W.Lin 
Page:  2831-2836 
Year  Published:  1997 
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Name  of  Conference:  28th  AIAA  Fluid  Dynamics  Conference/  4th  AIAA  Shear  Flow 
Control  Conference,  Snowmass  Village,  CO 

Title  of  Article:  “Example  of  output  regulation  for  distributed  parameter  systems” 
Author:  C.I.  B5rrnes,  D.S.  Gilliam,  A.  Isidori,  I.  Lauko,  and  V.I.  Shubov 
Year  Published:  1997. 


Books: 

Title  of  Book:  “Output  Regulation  of  Uncertain  Nonlinear  Systems” 

Author  of  Book:  C.I.Byrnes,  F.Delli  Priscoli,  A.  Isidori  Publisher:  Birkhauser  (Boston) 

Pages:  1-120 

Year  Published:  1997. 

Title  of  Book:  “Systems  and  Control  in  the  Twenty-First  Century” 

Editors  of  Book:  C.I.Byrnes,  B.N.  Datta,  D.S.  Gilliam,  C.F.  Martin 
Publisher:  Birkhauser  (Boston) 

Pages:  1-434 
Year  Published:  1997. 

Book  Chapters: 

Title  of  Book:  “Colloquim  on  Automatic  Control” 

Editor  of  Book:  C.Bonivento,  G.Marri,  G.Zanasi 

Title  of  Chapter:  “Semiglobal  robust  regulation  of  nonlinear  systems” , 

Author  of  Chapter:  A. Isidori 
Publisher:  Springer  Verlag 
Pages:  31-59,  Year  Published:  1996. 

Title  of  Book:  “Systems  and  Control  in  the  Twenty-First  Century” 

Editor  of  Book:  C.I.Byrnes,  B.N.  Datta,  D.S.  Gilliam,  C.F.  Martin 

Author  of  Chapter:  C.I.  Byrnes  and  A.  Lindquist 

Title  of  Chapter:  “On  duality  between  filtering  and  interpolation” 

Pages:  101-136  Publisher:  Birkhauser  (Boston) 

Title  of  Book:  “Current  and  Future  Directions  in  Applied  Mathematics” 

Editor  of  Book:  M.Alber,  B.  Hu,  J.  Rosenthal 

Author  of  Chapter:  C.I.  Byrnes,  H.J.  Landau,  and  A.  Lindquist 

Title  of  Chapter:  “On  the  well-posedness  of  the  rational  covariance  extension  problem” 

Pages:  83-108 

Publisher:  Birkhauser  (Boston) 
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5  Scientific  Interactions 

Our  prinicpal  investigating  team  has  enjoyed  collaborative  research  with  engineering  research 
and  development  personnel  at  the  McDonnell- Douglas  Aircraft  Division  of  Boeing  in  St. 
Louis,  MO,  at  Lockheed  -  Martin  in  Dallas,  TX  and  scientific  interaction  with  AFOSR 
personnel  at  Bolling  AFB,  Ft.  Eglin  AFB  and  Wright  Patterson  AFB.  Our  most  recent 
contacts  have  been; 

During  the  month  of  February,  1997,  Dr.  Christopher  I.  Byrnes  presented  a  technical 
briefing  on  nonlinear  control  systems  and  problems  in  flow  control  for  Dr.  Marc  Jacobs 
AFOSR  and  Dr.  Charles  Holland  AFOSR.  Other  participants  included  Dr.  Rowena  Eber- 
hardt  (McDonnell  Douglas  Aerospace,  St.  Louis,  MO),  Dr. David  E.  Paxekh  (McDonnell 
Douglas  Aerospace,  St.  Louis,  MO),  Dr.  Kevin  Wise  (McDonnell  Douglas  Aerospace,  St. 
Louis,  MO),  and  Dr.  David  S.  Gilliam,  Dr.  Ervin  Roden,  Dr.  Massoud  Amin,  Dr.  Barna 
Szabo,  Mr.  James  Ramsey  from  Washington  University. 

Professors  Christopher  I.  Byrnes  and  Alberto  Isidori  have  consulted  with  Dr.  Rowena  Eber- 
hardt  (McDonnell  Douglas  Aerospace,  St.  Louis,  MO)  on  robust  output  regulation  of  non¬ 
linear  systems.  This  work  is  aimed  at  extending  robust  output  regulation  methodologies  to 
fiight  controller  design  of  tailess  aircraft,  including  the  X36  and  MANX  aircraft. 

Alberto  Isidori  taught  a  course  on  Nonlinear  Systems  and  Control  at  the  NASA  Research 
Center  in  Langley,  VA.  Program.  Part  I  (15  hours):  nonlinear  systems  structure,  feedback 
linearization,  zero  dynamics.  Part  II  (15  hours):  nonlinear  system  stabilization,  regulation 
and  tracking.  Part  I  was  taught  in  November  1996,  Part  II  in  February  1997. 

During  the  month  of  April,  1997,  Dr.  David  Gilliam  consulted  with  researchers  in  the  Air 
Force  Center  for  Optimal  Design  and  Control  at  Virginia  Tech  &  State  University  concern¬ 
ing  the  design  of  feedback  control  laws  for  the  problem  of  output  regulation  of  distributed 
parameter  systems.  During  this  same  visit  Dr.  Gilliam  served  as  an  advisor  on  questions 
concerning  the  existence  of  global  dynamics  and  attractors  for  their  work  on  Thermal  Con¬ 
vective  Processes. 

We  have  presented  many  invited  lectures  and  colloquia  nationally  and  internationally: 
September  1994 

“Feedback  Stabilization  About  Attractors  and  Inertial  Manifolds,”  Invited  lecture  at 
Royal  Institute  of  Technology,  Stockholm,  lecture  presented  by  Professor  C.I.  Byrnes 

February  1995: 

“Feedback  Stabilization  About  Attractors  and  Inertial  Manifolds”  Universita  di  Roma 
-  La  Sapienza  Presented  by  Professor  C.  I.  Byrnes 

April,  1995: 

“Feedback  Stabilization  in  the  Large  for  Nonaffine  Nonlinear  Systems  with  Stable  Free 
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Dynamics,”  3rd  SIAM  Conference  on  Control  and  Its  Applications,  St.Louis,  Presented 
by  Professor  Wei  Lin. 

June  1995: 

“Nonlinear  Control  Systems,”  AFOSR  Contractors  Review  Meeting  Minneapolis-St. 
Paul,  Presented  by  Professor  C.I.  Byrnes 

“A  Riccati  Equation  for  Partial  Stochastic  Realization  Theory”  Royal  Institute  of 
Technology  Stockholm,  Sweden  Presented  by  Professor  C.I.  Byrnes 

“Some  New  Results  on  Stability  and  Observability  of  Discrete-time  Autonomous  Sys¬ 
tems”  American  Control  Conference  Seattle,  Presented  by  Professor  Wei  Lin. 

“Feedback  Stabilization  of  General  Nonlinear  Systems,”  IFAC  Nonlinear  Control  Sys¬ 
tems  Symposium,  Tahoe  City,  Presented  by  Professor  Wei  Lin. 

“On  the  Dynamics  of  Boundary  Controlled  Nonlinear  Distributed  Parameter  Systems” 
IFAC  Nonlinear  Control  Systems  Symposium,  Tahoe  City,  Presented  by  Professor  D. 
S.  Gilliam 

October  1995: 

“Robust  Regulation  of  Nonlinear  Systems,”  Invited  talk  presented  by  Prof.  C.I.  Byrnes 
at  VPI  &  State  University. 

“Robust  Regulation  of  Nonlinear  Systems,”  Plenary  talk  presented  by  Prof.  C.I.  Byrnes 
at  NCSU  Southeast  Regional  Differential  Equations  Conference. 

“On  the  Dynamics  of  Boundary  Controlled  Convective  Reaction  Diffusion  Equations,” 
Invited  talk  presented  by  D.S.  Gilliam  in  the  minisymposium  “Control  and  Compu- 
tatuional  Fluid  Dynamics”  at  the  1995  SIAM  Annual  Meeting  at  Charlotte,  North 
Carolina. 

“Electromagnetic  Waves  Over  a  Three  Dimensional  Conducting  Half  Space  Invited 
lecture  presented  by  D.S.  Gilliam  in  the  Center  For  Research  in  Scientific  Computa¬ 
tion  at  North  Carolina  State  University. 

“Low-Gain  Feedback  Stabilization  of  Controllable  Nonlinear  Systems  with  Stable  Free 
Dynamics”  Invited  lecture  presented  by  Prof.  Wei  Lin  at  the  Center  for  Control  Engi¬ 
neering  and  Computation,  Dept,  of  Electrical  h  Computer  Engineering,  University  of 
California,  Santa  Barbara. 

November  1995: 

“Robust  Regulation  of  Nonlinear  Systems,”  Invited  talk  presented  by  Prof.  C.I.  Byrnes 
in  the  Departemnt  of  Electrical  Engineering  at  the  University  of  Hong  Kong,  Hong 
Kong. 
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“Output  Regulation  Revisited”  Distinguished  lectures  series,  Department  of  Electrical 
Engineering,  IOWA  State  University,  presented  by  Dr.  Alberto  Isidori. 

“Robust  Stabilization  and  Adaptive  Regulation  of  Minimum-Phase  Nonlinear  Systems 
with  Uncertainty”  Invited  lecture  presented  by  Prof.  Wei  Lin  at  the  Dept,  of  Electrical 
Engineering,  Iowa  State  University,  Ames. 

December  1995: 

“Disturbance  Attenuation  for  a  Class  of  Nonlinear  Non-Minumum  Phase  Systems,” 
Invited  lecture  presented  by  Prof.  Alberto  Isidori  at  the  34th  Control  and  Decision 
Conference. 

“High  Gain  Limit  for  Boundary  Controlled  Convective  Reaction  Diffusion  Equations” 
lecture  presented  by  D.S.  Gilliam  at  the  34th  Control  and  Decision  Conference. 

“Well-Posedness  for  Controlled  Nonlinear  Damped  Membranes  with  Fixed  Boundary” 
lecture  presented  by  D.S.  Gilliam  at  the  34th  Control  and  Decision  Conference. 

“Global  Robust  Stabilization  of  Minimum-Phase  Nonlinear  Systems”  lecture  presented 
by  W.Lin  at  the  34th  Control  and  Decision  Conference. 

“Global  Stabilization  of  Discrete  Non-Affine  Systems”  lecture  presented  by  W.Lin  at 
the  34th  Control  and  Decision  Conference. 

“Mixed  H2/H- infinity  Control  of  Nonlinear  Systems”  Author:  W.  Lin  lecture  presented 
by  W.Lin  at  the  34th  Control  and  Decision  Conference. 

January  1996: 

“Dynamics  of  systems  governed  by  convection  reaction  diffusion  equations,”  lectmre 
presented  at  the  Joint  AMS/MAA  Annual  Meeting  in  Orlando,  Florida. 

“The  Regulator  Problem  for  Linear  Distributed  Parameter  Systems,”  Invited  lecture 
presented  by  D.S.  Gilliam  in  the  Department  of  Systems  Science  and  Math  at  Wash¬ 
ington  University  in  St.  Louis. 

March  1996: 

“Complete  Parameterization  of  Solutions  to  the  Rational  Covariance  Problem,”  Invited 
lecture  presented  by  C.I.  Byrnes  at  Ecole  Normale  Superieure,  Paris. 

April  1996: 

“Complete  Parameterization  of  Solutions  to  the  Rational  Covariance  Problem,”  Ple¬ 
nary  lecture  presented  by  C.I.  Byrnes  at  Conference  on  the  Future  Directions  in  Applied 
Mathematics  at  Notre  Dame. 
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June  1996: 

“Global  Stabilization  of  Discrete-Time  Nonlinear  Systems  via  Bounded  Feedback”  In¬ 
vited  talk  presented  by  Prof.  Wei  Lin  at  the  96  Mathematical  Theory  of  Networks  and 
Systems  (MTNS)  Conference,  St.  Louis,  MO. 

“The  Structure  of  Attractors  for  a  Boundary  Controlled  Viscous  Burgers’  Equation,” 
Invited  talk  presented  by  Andras  Balogh  at  the  96  Mathematical  Theory  of  Networks 
and  Systems  (MTNS)  Conference,  St.  Louis,  MO. 

“Recent  Results  on  Boundary  Control  for  a  Class  of  Higher  Dimensional  Convection- 
Reaction  Diffusion  Equations,”  Invited  talk  presented  by  David  S.  Gilliam  at  the  96 
Mathematical  Theory  of  Networks  and  Systems  (MTNS)  Conference,  St.  Louis,  MO. 

’’Nonlinear  Observers  for  Hopf  Bifurcating  Systems”  Invited  talk  presented  by  Dr. 
V.  Sundarapandian  at  the  96  Mathematical  Theory  of  Networks  and  Systems  (MTNS) 
Conference,  St.  Louis,  MO. 

“Adaptive  Control  of  Minimum-Phase  Nonlinear  Systems”  presented  by  Prof.  Wei 
Lin  at  the  96  Mathematical  Theory  of  Networks  and  Systems  (MTNS)  Conference,  St. 
Louis,  MO. 

“Robust  Nonlinear  Output  Regulation,”  Invited  talk  presented  by  Prof.  Alberto  Isidori 
at  the  96  Mathematical  Theory  of  Networks  and  Systems  (MTNS)  Conference,  St. 
Louis,  MO. 

“Passivity,  Bounded  State  Feedback  and  Global  Stabilization  of  Nonlinear  Systems” 
presented  by  Prof.  Wei  Lin  at  the  13th  IFAC  World  Congress. 

“Solutions  to  the  Output  Regulation  Problem  of  Linear  Singular  Systems,”  presented 
by  Prof.  Wei  Lin  at  the  13th  IFAC  World  Congress. 

August  1996: 

“On  the  Rational  Caratheodory  Extension  Problem,”  Invited  lextrue  presented  by 
Prof.  C.I.  Byrnes  at  the  Fifth  Bozeman  Conference  on  Computation  and  Control. 

“The  Regulator  Problem  for  Linear  Distributed  Parameter  Systems,”  Invited  lecture 
presented  by  Istvan  Lauko  at  the  Fifth  Bozeman  Conference  on  Computation  and  Con¬ 
trol. 

“On  the  Dynamics  and  Control  of  a  Class  of  Higher  Dimensional  Convection  Reac¬ 
tion  Diffusion  Equations,”  Invited  lecture  presented  by  David  Gilliam  at  the  Fifth 
Bozeman  Conference  on  Computation  and  Control. 
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“Numerical  Stationary  Solutions  for  Burgers’  Equation  with  Neumann  Boundary  Con¬ 
ditions,”  Invited  lecture  presented  by  David  Gilliam  at  the  Fifth  Bozeman  Conference 
on  Computation  and  Control. 

“Regularization  and  Control  of  Three  Dimensional  Navier  Stokes  Equations,”  Invited 
lecture  presented  by  Andras  Balogh  at  the  Fifth  Bozeman  Conference  on  Computation 
and  Control. 

October  1996: 

“Zeros  and  Zero  Dynamics  for  Linear  Distributed  Parameter  Systems,”  A  series  of  lec¬ 
tures  given  by  Dr.  David  Gilliam  in  the  Texas  Tech  University  Joint  Civil  Engineering, 
Mathematics  and  Mechanical  Engineering  Seminar. 

December  1996: 

“Global  normal  forms  for  MIMO  nonlinear  systems,  with  application  to  stabilization 
and  disturbance  attenuation”  Lecture  presented  by  Dr.  Alberto  Isidori  at  the  35rd 
IEEE  Conf.  Decision  and  Control  (Kobe,  Japan). 

“L2  disturbance  attenuation  and  performance  bounds  for  linear  non-minimum  phase 
square  invertible  systems”  Lecture  presented  by  Dr.  Alberto  Isidori  at  the  35rd  IEEE 
Conf.  Decision  and  Control  (Kobe,  Japan). 

“Immersions  and  the  internal  model  principle:  tools  for  robust  nonlinear  control”  Lec¬ 
ture  presented  by  Dr.  Christopher  1.  Byrnes  at  the  35rd  IEEE  Conf.  Decision  and 
Control  (Kobe,  Japan). 

“The  rational  covariance  extension  problem  with  applications  to  speech  synthesis” 
Invited  lecture  presented  by  Dr.  Christopher  1.  Byrnes  at  The  University  of  Tokyo. 

“Robust  nonhnear  control”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes 
at  Dewan  Riset  Nasional,  Serpong,  Indonesia. 

January  1997: 

“Global  L2-Gain  State  Feedback  Design  for  a  Class  of  Nonlinear  Systems”  Lecture 
presented  by  Dr.  Alberto  Isidori  at  COSY  Workshop  on  Nonlinear  Systems  held  at 
the  ETH-Zurich,  Switzerland. 

February  1997: 

“The  geometry  of  positive  real  functions  with  applications  to  signal  processing  and  to 
speech  synthesis”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at  The  Texas 
Tech  University. 

March  1997: 

“The  rational  covariance  extension  problem  with  applications  to  speech  synthesis”  In¬ 
vited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at  Texas  Tech  University. 
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“Robust  Output  Regulation  of  Nonlinear  Systems”  Invited  lecture  presented  by  Dr. 
Alberto  Isidori  at  the  University  of  Aalborg,  Denmark. 

April  1997: 

“On  the  Global  Dynamics  of  a  Class  Boundary  Controlled  Hydrodynamic  Systems  in 
Higher  Dimensions,”  Lecture  presented  by  Dr.  Victor  I.  Shubov  at  the  Texas  Partial 
Differential  Equations  Conference  in  Denton,  TX. 

“Harmonic  forcing  for  linear  distributed  parameter  systems”  Invited  lecture  presented 
by  Istvan  Lauko  in  the  Texas  Tech  University  Joint  Civil  Engineering,  Mathematics 
and  Mechanical  Engineering  Seminar. 

May  1997: 

“Nonlinear  control  systems”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at 
the  1997  AFOSR  Contractors  Review  Meeting,  Wright  Patterson  AFB,  Dayton,  OH. 

June  1997: 

“On  the  covariance  extension  problem”  5th  IEEE  Mediterranian  Conference,  Paphos, 
Cyprus  Invited  lecture  presented  by  Professor  C.  I.  Byrnes 

“Robust  Output  Regulation  of  Nonlinear  Systems”  Invited  lecture  presented  by  Dr. 
Alberto  Isidori  at  the  Imperial  College,  London,  UK. 

“Robust  Nonlinear  Control”  Invited  lecture  presented  by  Dr.  Christopher  1.  Byrnes  at 
the  28th  AIAA  Fluid  Dynamics  Conference/  4th  AIAA  Shear  Flow  Control  Conference, 
Snowmass  Village,  CO 

July  1997: 

Dr.  David  S.  Gilliam  participated  in  12th  International  Conference  of  the  American 
Meterological  Society  on  Boundary  Layers  and  Turbulence,  Vancouver,  BC. 

“On  Duality  between  Filtering  and  Interpolation”  University  of  Tokyo,  Japan  Invited 
lecture  presented  by  Professor  Anders  Lindquist. 

December  1997: 

“Global  L2-gain  State  Feedback  Design  for  a  Class  of  Nonlinear  Systems,”  36th  IEEE 
Conf.  Decision  and  Control  (San  Diego,  CA,  December  1997),  Lecture  presented  by 
Professor  A. Isidori 

“Harmonic  Forcing  for  a  Class  of  Nonlinear  Systems,”  36th  IEEE  Conf.  Decision  and 
Control  (San  Diego,  CA,  December  1997),  Lecture  presented  by  Professor  A. Isidori 

“Computational  Methods  for  Feedback  Control  of  Distributed  Parameter  Systems,” 
Special  Session  at  36th  IEEE  Conf.  Decision  and  Control  (San  Diego,  CA,  December 
1997),  Organizer  and  Chair  Professor  David  Gilliam. 
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New  Discoveries 


One  of  the  highlights  of  our  research  program  during  this  reporting  period  has  been  the 
research,  the  preparation  and  the  publication  of  our  recent  monograph  describing  a  solution 
to  the  problem  of  robust  output  regulation  for  lumped  nonlinear  systems.  Since  systems 
containing  unknown  real  parameters  are  rarely  detectable  or  observable,  when  the  state  is 
augmented  to  include  the  unknown  parameters,  it  is  typically  impossible  to  directly  com¬ 
pute  a  (for  example,  stablizing)  feedback  law  as  a  function  of  the  system  output.  A  key 
development  has  been  the  identification  of  an  alternative  to  state  feedback/state  observer 
design  which  is  better  suited  for  robust  control.  This  alternative  is  based  on  a  combination 
of  the  internal  model  principle  and  the  notion  of  immersion.  Using  these  tools,  a  solution  of 
the  robust  output  regulation  problem  has  been  obtained.  One  early  example  of  an  alternative 
to  directly  computing  a  gain  from  more  data  than  is  actually  need  to  produce  the  gain  is 
the  so-called  fast  filtering  algorithm  for  Kalman  filtering.  Thus,  instead  of  computing  the 
Kalman  gain  as  the  product  of  a  positive  solution  of  a  Riccati  equation  and  a  state  estimate, 
fast  filtering  propagates  the  Kalman  gain  as  a  vector,  along  with  a  covector. 

Another  especially  important,  but  unanticipated,  success  in  our  research  effort  has  begun 
with  the  development  of  the  phase  portrait  for  the  fast  filtering  algorithm  for  scalar  systems, 
as  a  nonlinear  dynamical  system.  This  phase  portrait,  and  its  relation  to  the  geometry 
of  spaces  of  positive  real  functions,  has  recently  been  used  to  solve  a  long  standing  open 
problem  in  signal  processing  and  speech  synthesis,  with  important  implications  in  the  digital 
encoding  and  transmission  of  human  speech.  Indeed,  we  have  developed  a  geometric  duality 
between  filtering  and  interpolation,  tying  together  two  fundamental  classes  of  problems, 
rational  covariance  extension  and  the  dynamics  of  fast  filtering  algorithms.  Perhaps  the  most 
important  accomplishment  during  this  particular  research  period  has  been  the  discovery  of 
a  convex  minimization  problem  that  leads  to  an  effective  algorithm  for  speech  synthesis  and 
speech  analysis.  These  algorithms,  when  coded  on  a  programmable  digital  signal  processing 
chip  and  coupled  with  A/D  and  D/A  converters,  allow  for  the  encoding  and  decoding, 
and  hence  the  transmission,  of  human  speech.  In  addition  to  wireless  speech  transmission 
technologies,  there  are  also  potential  applications  to  secure  communications  systems.  A 
patent  for  using  this  methodology  in  speech  synthesis  applications  is  pending. 

Among  our  most  significant  accomplishments  in  our  research  in  the  area  of  distributed 
parameter  systems  is  the  development  of  geometric  theory  of  output  regulation  for  linear 
distributed  parameter  systems  with  bounded  actuators  and  sensors.  In  particular,  for  a  very 
general  class  of  distributed  parameter  plants  with  bounded  input  and  output  operators  and 
for  reference  signals  and  disturbances  generated  by  finite  dimensional  exogenous  systems,  we 
have  given  necessary  and  sufficient  conditions  for  solvability  of  the  state  feedback  regulator 
problem  (SFRP)  and  error  feedback  regulator  problem  (EFRP).  These  conditions  are  given 
in  terms  of  the  solvability  of  a  pair  of  regulator  equations  which  are  easily  solved  numerically 
off-line.  The  solution  of  these  equations  provide  explicit  and/or  numerical  solutions  which 
yield  an  exact  and/or  approximate  feedback  law  for  the  SFRP  or  controller  for  solving  the 
EFRP.  These  techniques  have  been  successfully  applied  to  problems  in  asymptotic  tracking  of 
harmonic  reference  signals  and  set  point  control  for  both  parabolic  and  hyperbolic  examples. 
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Further,  for  the  problem  of  tracking  a  periodic  motion  and  for  plants  whose  state  dynamics 
are  governed  discrete  Riesz  spectral  operators  we  have  been  able  to  characterize  solvability 
of  the  regulator  equation  in  terms  of  certain  nonreasonance  conditions  related  to  the  poles 
and  zeros  of  the  open  loop  transfer  function.  For  systems  for  which  the  zero  dynamics 
exists,  we  are  also  able  to  characterize  solvability  of  the  regulator  equations  in  terms  of  the 
existence  of  an  isomorphic  copy  of  the  exosystem  in  the  composite  systems  zero  dynamics, 
complimentary  to  the  zero  dynamics  of  the  plant.  Sponsored  work  with  McDonnell  Douglas 
Aircraft  is  underway  to  transition  output  regulator  technologies  to  the  problem  of  shaping 
the  airflow  downstream  from  a  wing  section  in  order  to  either  produce  more  lift  for  tail  flaps, 
or  to  better  mix  the  jet  exhaust  in  an  effort  to  lower  heat  and  infrared  signals. 

7  Additional  Information 

Our  research  team  has  been  involved  in  the  following  transitions  of  technology.  Dr.  Christo¬ 
pher  I.  B3n:nes  and  Dr.  David  S.  Gilliam  together  with  Dr.  David  Parekh  (McDonnell 
Douglas  Aerospace,  St.  Louis,  MO)  have  made  progress  on  robust  nonlinear  control  design 
for  problems  of  flow  control.  The  intended  application  being  the  development  of  a  robust 
control  strategy  for  the  control  of  airflow  to  increase  lift  on  tail  flaps  on  aircraft  such  as  the 
C17.  In  May  1997,  Drs.  Christopher  I.  Byrnes  and  Anders  G.  Lindquist  have  submitted, 
through  Washington  University,  the  patent  discloure,  “A  new  method  and  device  for  speech 
synthesis,”  to  the  U.S.  Patent  Office. 
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